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ABSTRACT 


Objective  Analysis  of  Tropical  Cyclone  Intensity,  Strength,  and  Size 
Using  Routine  Aircraft  Reconnaissance  Data.  (May  1986) 

Charles  Baynard  Stanfield,  B.S.,  Florida  State  University 
Chairman  of  Advisory  Committee:  Dr.  Aylmer  H.  Thompson 

\he  feasibility  of  objectively  analyzing  routine  aircraft  recon¬ 
naissance  data  for  the  purpose  of  quantifying  tropical  cyclone  inten¬ 
sity,  strength,  and  size  is  examined.  A  computer  program  is  developed 
which  may  be  used  in  near  real  time  or  after  the  fact  to  evaluate 
localized  pressure/wind  relationships  in  the  tropical  cyclone  environ¬ 
ment.  This  program  compensates  for  the  system  motion  and  the  relative 
position  of  the  point  of  observation  relative  to  the  vortex  center  loca 
tion  at  flight  level  and  at  the  surface  (thus  accounting  for  the  verti¬ 
cal  tilt  of  the  center).  A  representative  set  of  data  is  obtained  over 
a  13  month  period  for  the  entire  spectrum  of  storms  from  tropical  de¬ 
pression  to  super  typhoon.  These  data  are  used  to  try  to  establish 
empirical  pressure/wind  relationships  and  a  means  of  determining  effec¬ 
tive  storm  size.  It  is  shown  that  a  program  of  this  nature  may  be  used 
with  gradient  wind  and  pressure  gradient  relationships  to  evaluate  in¬ 
tensity  and  strength  and  to  define  storm  size,  provided  adequate  data 
are  available  at  sufficient  distances  from  the  center. 
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CHAPTER  I 


INTRODUCTION 

During  the  early  years  of  tropical  cyclone  reconnaissance,  the  lack 
of  technology  limited  the  ability  of  observers  aboard  reconnaissance 
aircraft  to  collect  meaningful  information  in  the  periphery  of  the 
tropical  cyclone.  The  rather  undisciplined  approach  of  reconnaissance 
procedures  led  scientists  using  the  data  to  use  only  a  few  variables 
(i.e.,  minimum  sea  level  pressure)  in  developing  their  equations  and 
models  of  tropical  cyclone  structure. 

With  advances  in  technology,  the  quality  and  density  of  data  in  the 
tropical  cyclone  environment  increased.  The  content,  format,  and 
collection  of  routine  tropical  cyclone  reconnaissance  information  have 
become  more  systematic  and  disciplined  in  recent  years.  Still,  except 
for  data  collected  by  the  United  States  National  Oceanic  and  Atmospheric 
Administration  Research  Flight  Facility,  little  use  is  made  of  peri¬ 
pheral  data  that  are  being  collected  on  routine  reconnaissance  missions. 
With  the  use  of  these  data,  the  structure  of  specific  tropical  cyclones 
in  various  stages  of  development  can  be  understood  better.  In  addition 
to  the  intensity  (maximum  sustained  wind  of  the  storm),  the  effects  of 
overall  strength  (average  wind  speed)  and  size  (radius  of  the  circula¬ 
tion)  of  a  specific  storm  can  be  determined  more  accurately  in  real  time 
if  concrete  relationships  with  pressure  gradient  are  established. 

The  citations  on  this  and  the  following  pages  follow  the  style  of 
the  Monthly  Weather  Review. 
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Estimates  of  surface  winds  by  observers  in  aircraft  flying  at  the 
70  kPa  pressure  level  (approximately  3000  m  above  sea  level)  and  even  at 
low  level  (less  than  500  m  above  sea  level),  have  been  neglected  as  re¬ 
liable  depictions  of  the  actual  wind  profile  of  the  tropical  cyclone. 
Reasons  for  this  include:  1)  individual  observer  misinterpretation,  2) 
observer  inexperience,  3)  observations  taken  at  non-specified  levels  for 
surface  winds,  4)  observations  taken  in  poor  visibility  conditions,  5) 
the  relatively  broad  classifications  of  surface  wind  speed,  and  6)  the 
inability  to  see  the  entire  sea  surface  area  of  the  tropical  cyclone. 
This  last  reason  includes  the  factors  of  intervening  clouds  and  restric¬ 
ted  flight  in  the  area  of  the  storm.  Dynamically,  the  maximum  winds  are 
associated  with  the  eye  wall  cloud  and  convective  feeder  bands.  In 
these  regions,  the  sea  surface  usually  is  obscured  by  cloud,  and  large 
amounts  of  ocean  water  are  picked  up  and  carried  in  the  air.  Because  of 
these  obstacles,  maximum  surface  wind  may  be  estimated  incorrectly  by 
visual  observation.  Because  Doppler  radar  may  track  airborne  moisture 
as  well  as  the  sea  surface,  the  flight  level  wind  may  be  estimated  in¬ 
correctly  by  instruments. 

Even  though  problems  still  exist  with  surface  wind  estimates  by 
aircraft  observers,  more  accurate  definitions  of  wind  speeds  below  67 
m  s  ^  have  resulted  because  of  the  years  of  experience  since  the  first 
reconnaissance  missions  were  flown.  This  study  will  bring  out  the  weak¬ 
nesses  and  strengths  of  using  these  data,  and  recommendations  will  be 
offered  as  to  how  to  improve  the  usefulness  of  the  data. 


CHAPTER  II 


OBJECTIVE 

This  research  examines  the  potential  of  routine  aircraft  reconnais¬ 
sance  observations  to  determine  objectively  tropical  cyclone  intensity, 
strength,  and  size  characteristics.  The  term  "routine"  refers  to  those 
means  of  data  collection  used  on  U.S.  Air  Force  aircraft  reconnaissance 
missions  of  the  late  1970's  and  early  1980‘s. 

The  primary  objective  is  to  derive  empirical  relationships  (if 
possible)  that  relate  the  observed  70  kPa  height  gradient  or  the  surface 
pressure  gradient  to  the  intensity,  strength,  and  size  characteristics 
of  tropical  cyclones  at  the  earth's  surface.  First,  relationship  of  the 
70  kPa  height  gradient  or  sea  level  pressure  gradient  to  the  observed 
winds  is  determined  using  a  relatively  coarse  distribution  of  aircraft 
observations  centered  in  relation  to  the  vortex  center.  Second,  the 
gradient  wind  is  computed  using  the  height/pressure  gradients  and  com¬ 
pared  to  the  observed  winds.  The  effects  of  system  motion  and  vertical 
tilt  of  the  vortex  center  are  accounted  for  in  the  program.  A  determi¬ 
nation  is  made  of  which  method  (the  empirical  relationships  of  observed 
winds  and  pressure  gradients  or  simply  the  gradient  wind  relationships) 
is  best  in  defining  the  location  and  magnitude  of  the  maximum  wind  zone 
(intensity)  and  the  overall  wind  speed  profile  (strength)  of  tropical 
cyclones.  In  addition,  a  method  is  examined  for  defining  storm  size 
using  an  assumed  pressure  profile  in  data  sparse  areas. 

The  results  are  interpreted  in  light  of  the  data  gathering  techni¬ 
ques  of  the  past  as  well  as  the  future,  and  recommendations  are  made  for 
upgrading  of  the  aircraft  reconnaissance  platform. 
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CHAPTER  III 
BACKGROUND 

The  endeavor  to  establish  a  relationship  between  tropical  cyclone 
winds  and  pressure  began  in  1939.  Using  data  from  Pacific  typhoons, 
Takahashi  (1939)  was  the  first  to  develop  the  rather  fundamental  rela¬ 
tionship, 

=  '3.4  (p„  -  .  (1) 

from  the  cyclostrophic  wind  equation.  This  equates  the  maximum  surface 

wind  speed  in  knots  (V^)  to  the  difference  between  the  cyclone's  minimum 

sea  level  pressure  (p  )  and  the  environmental  pressure  (p  )  in  milli- 

ri 

bars.  Most  research  since  Takahashi  has  been  only  to  modify  this  rela¬ 
tionship. 

McKnown  (1952)  and  his  assistants  used  aircraft  reconnaissance  data 
from  Pacific  typhoons  and  further  modified  the  equation.  His  equation 
is 

=  (20  -  «/5)  (1010  -  p^)°'5  ,  (2) 


where  4>  is  latitude  in  degrees. 

Fletcher  (1955)  used  data  collected  by  the  Corps  of  Engineers  in 
the  hurricane  of  26-27  August,  1949  at  Lake  Okeechobee,  Florida  to  re¬ 
vise  Takahashi 's  equation,  obtaining 


V_  =  16  (p_  -  p.) 
m  '^n  '^c' 


•  •  -  •  •  «  *  . 


V  "W  s'  s*  s  s 


Here  was  defined  as  the  environmental  pressure  at  the  "outer  edge"  of 
the  cyclone  and  was  assumed  to  be  either  1010  mb  or  2-3  mb  higher  than 
the  outermost  closed  surface  isobar,  as  analyzed  on  the  surface  synoptic 
chart. 

Prior  to  1956,  equation  development  was  based  upon  the  peak  gusts 
of  the  storm  rather  than  the  sustained  (one  minute  average)  surface  wind 
in  knots  (V^).  Kraft  (1961)  was  the  first  to  modify  Fletcher's  equation 
to  account  for  this.  He  used  a  slightly  different  value  for  the  Atlan¬ 
tic  environmental  pressure  to  get 

=  14  (1013  -  p^)°'^  .  (4) 

This  equation  is  used  now  in  the  Atlantic  area. 

With  the  increase  of  reconnaissance  flights  at  70  kPa,  an  extra¬ 
polation  technique  was  developed  to  estimate  the  central  sea  level 
pressure  of  tropical  cyclones  from  70  kPa  data  (Jordan,  1957).  Modifi¬ 
cations  were  made  to  compute  maximum  winds  based  upon  minimum  70  kPa 
heights  alone.  These  equations  were  used  throughout  the  entire  1960's 
and  most  of  the  1970's. 

Atkinson  and  Holliday  (1977)  were  the  first  to  examine  surface  wind 
reports  from  exposed  coastal  or  island  stations  where  there  was  a  high 
probability  that  the  station  actually  experienced  the  maximum  winds  of 
the  storm.  They  developed  an  empirical  relationship  that  resembled  the 
cyclostrophic  equation,  yet  also  represented  the  best  nonlinear  fit  to 
the  data.  Their  equation  may  be  written  as 

0.6444 


V  .  =  6.7  (1010  -  pj 


(5) 


Even  though  (5)  has  proven  to  be  a  very  good  equation  when  applied 
to  typhoons  of  the  western  North  Pacific,  there  are  some  unanswered 
problems.  Exact  sustained  one-minute-average  wind  speeds  are  very  dif¬ 
ficult  to  obtain  from  wind  recorder  data  because  of  the  rapid  variation 
of  wind  speed  during  high  speed  conditions.  Peak  gusts  were  used  and 
adjusted  to  sustained  wind  speed  values  using  adjustment  factors  de¬ 
veloped  by  the  Air  Force  Cambridge  Research  Laboratory  (Sissenwine  e^ 
al . ,  1973).  However,  these  adjustment  factors  were  not  developed  for 
tropical  cyclone  wind  and  turbulence  cases.  In  addition,  (5)  is  the 
representation  of  the  nonlinear  regression  line  which  best  fits  the  data 
and,  therefore,  represents  a  mean  of  the  data.  Observations  above  the 
regression  line,  including  those  from  sensors  that  broke  or  blew  away  at 
or  before  the  maximum  winds  were  recorded,  are  indicative  of  higher 
winds  due  to  local  pressure  gradient  or  supergradient  relationships. 

Analytic  models  were  developed  by  Depperman  (1947),  Schloemer 
(1954),  and  Holland  (1980)  to  reproduce  the  horizontal  profiles  of  tro¬ 
pical  cyclone  pressure  and  wind.  The  models  contain  parameters  which 
are  estimated  empirically  from  climatological  composites  or  from  obser¬ 
vations  of  specific  storms.  From  these  parameters,  a  geometric  rela¬ 
tionship  is  derived  for  winds  both  inside  and  outside  the  radius  of 
maximum  wind.  Although  analytic  models  are  useful  in  approximating  the 
general  wind  profile  of  the  average  hurricane/typhoon,  the  wind  and 
pressure  profiles  of  tropical  cyclones  vary  considerably  (Gray,  1981). 

As  mentioned  earlier,  a  tropical  cyclone  can  be  analyzed  in  three 
ways  (Merrill,  1982).  The  "intensity"  is  defined  as  the  maximum  sus¬ 
tained  wind  speed  of  the  storm,  and  even  though  it  can  realistically  be 


approximated  for  stronger  storms  from  minimum  sea  level  pressure  alone, 
the  location  and  extent  of  this  maximum  wind  band  can  only  be  evaluated 
with  sufficient  data.  Storm  "strength"  is  defined  as  the  average  wind 
speed  observed  within  the  radial  limits  of  the  storm,  and  storm  "size" 
is  defined  as  a  measurement  of  this  radius.  Little  attention  has  been 
given  to  these  latter  two  aspects.  In  addition,  most  research  and 
modeling  efforts  have  been  concentrated  on  tropical  cyclones  which  have 
attained  hurricane/ typhoon  strength  of  33  m  s”^ .  Comparatively  little 
is  known  about  the  highly  variable  tropical  depression  and  tropical 
storm  structure,  even  though  most  tropical  cyclones  either  never  attain 
hurricane/ typhoon  strength  or  remain  at  "depression"  or  "storm"  inten¬ 
sity  for  a  significant  portion  of  their  life  cycle. 

The  strength  of  a  storm  varies  considerably  with  time.  The 
strength  is  considered  to  be  dependent  on  the  overall  average  wind  speed 
within  the  cyclone,  not  just  the  maximum  wind.  To  demonstrate  this 
variability,  we  can  compare  the  strengths  of  two  super  typhoons  (Judy 
and  Tip)  of  the  western  North  Pacific  in  1979.  Super  typhoons  are  a 
subset  of  typhoons  in  which  the  maximum  surface  winds  reach  or  exceed 
67  m  s"^  at  some  time  during  the  life  of  the  storm.  Aircraft  reconnais¬ 
sance  data  collected  on  these  two  storms  indicated  that  Judy  had  a  much 
more  compact  pressure  and  wind  distribution  than  Tip.  This  comparison 
can  only  be  estimated  qualitatively,  since  adequate  quantitative  data 
during  the  periods  of  maximum  intensity  of  these  two  storms  could  not  be 
obtained  at  close  range  due  to  the  storm's  remote  location,  size,  and 
avoidance  by  all  but  reconnaissance  aircraft.  From  analysis  of  the 
reconnaissance  data,  however,  it  appears  that  Judy  had  an  extremely 


compact,  maximum  wind  band  with  all  the  significant  winds  within  a  very 
small  radius,  thus  affecting  only  a  small  area,  though  the  maximum  wind 
speed  was  higher  for  Judy  than  for  Tip.  However,  Tip  had  a  very  broad 
band  of  significant  winds,  which  affected  a  very  large  area  for  a  con¬ 
siderable  period  of  time.  Although  minimum  sea  level  pressure  and  maxi¬ 
mum  surface  winds  have  a  statistical  correlation,  the  wind  field  depends 
on  the  pressure  gradient,  not  the  minimum  pressure  itself.  An  objective 
method  to  depict  accurately  the  pressure  gradient  and  wind  profile  of  a 
tropical  cyclone  would  be  of  great  help  to  forecasters,  especially  in 
sparse  data  regions. 

In  addition  to  strength,  the  effects  of  the  cyclone's  size  are  con¬ 
siderable.  This  is  true  even  for  the  range  of  sizes  that  has  been  ob¬ 
served  for  tropical  cyclones  in  the  same  intensity  category.  Super 
Typhoon  Judy  attained  a  minimum  sea  level  pressure  of  88.7  kPa  and  had  a 
diameter  (d)  of  no  greater  than  1150  km,  as  determined  by  analysis  of 
its  outermost  closed  surface  isobar.  The  method  of  computing  the 
average  radius  of  the  outermost  closed  isobar  (AROCI)  is  shown  in  Fig. 

1.  According  to  Dunnavan  and  Diercks  (1980),  Super  Typhoon  Tip  reached 
a  record  minimum  sea  level  pressure  of  87  kPa  (only  slightly  lower  than 
that  of  Judy)  and  had  a  record  diameter  of  nearly  2220  km,  about  twice 
the  corresponding  value  for  Judy.  Even  though  a  more  extreme  example  of 
size  difference  could  be  cited  by  using  Hurricane  Tracy  of  Australia  in¬ 
stead  of  Typhoon  Judy,  the  comparison  of  Judy  and  Tip  (the  same  year  and 
general  location)  shows  the  vast  difference  in  size  of  certain  storms 
with  nearly  the  same  minimum  pressure  and  estimated  intensity.  The  area 
affected  by  cyclonic  winds  (0.25Trd  )  was  nearly  four  times  greater  for 


AROCI  =  (r^+r2+r2+r^)/4 


1.  Method  of  computing  the  average  radius  of  the 
outermost  closed  isobar  (AROCI).  On  most  maps 
the  isobars  are  at  a  4  mb  interval. 


Tip  than  for  Judy. 
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Bates  (1977)  developed  a  normalized  profile  of  vertical  wind  speed 
variation  for  the  standard  tropical  cyclone.  This  profile  provides  a 
theoretical  relationship  between  the  winds  at  the  70  kPa  level  and  those 
at  the  surface.  Consequently,  comparisons  of  winds  at  a  70  kPa  flight 
level  with  those  at  the  surface  would  be  expected  to  yield  fairly  con¬ 
sistent  results.  Whenever  surface  wind  speed  estimates  are  in  excess  or 
well  below  these  expectations,  the  estimates  are  usually  disregarded. 
However,  not  only  may  the  measured  flight  level  winds  be  in  error,  but 
the  normalized  wind  profile  is  meaningful  only  in  the  event  that  the 
storm  is  symmetric  and  vertically  aligned.  Huntley  and  Diercks  (1981) 
found  that,  in  1979,  47%  of  the  western  North  Pacific  typhoons  exhibited 
a  significant  vertical  tilt  between  their  surface  and  70  kPa  centers. 

The  concept  of  a  tilted  vortex  has  been  used  in  the  two-layer  discrete 
vortex  model  of  Khandekar  and  Rao  (1971)  to  study  the  short-term  dis¬ 
placements  of  tropical  vorticies  due  to  mutual  interaction.  During  the 
developing  stage  of  a  tropical  cyclone,  separation  between  70  kPa  and 
surface  centers  has  been  observed  to  be  over  100  km  and  on  rare  occa¬ 
sions  over  200  km.  In  the  later  stages,  the  tropical  cyclone  often 
tilts  with  height  significantly  in  the  direction  (downstream)  of  strong 
upper  level  steering  flow.  During  these  periods,  surface  and  upper 
level  winds  may  differ  considerably. 

Recent  studies  conducted  at  Colorado  State  University  by  Weather¬ 
ford  and  Gray  (1984),  Weatherford  (1985),  Merrill  (1985),  and  Edson 
(1985)  indicate  that  changes  in  tropical  cyclone  intensity  are  not 
correlated  well  with  changes  in  the  outer  wind  strength  (or  that  wind 


n 

between  about  100  to  280  km  from  the  center),  A  somewhat  better  corre¬ 
lation  exists  when  both  eye  size  and  minimum  pressure  are  known,  but  eye 
size  cannot  always  be  determined.  For  this  reason,  as  well  as  the  other 
factors  previously  discussed,  other  methods  of  data  analysis  are  de¬ 


sirable. 


CHAPTER  IV 


DATA  DESCRIPTION 

The  54th  Weather  Reconnaissance  Squadron  of  the  United  States  Air 
Force  is  the  primary  aircraft  reconnaissance  unit  to  gather  information 
from  the  numerous  tropical  cyclones  of  the  western  North  Pacific.  In  an 
average  year,  approximately  30  tropical  cyclones  develop  within  the 
western  North  Pacific  basin.  The  author  participated  in  reconnaissance 
flights  in  this  area  during  the  period  1979-1981. 

Over  one  full  year  of  reconnaissance  data  was  used  in  this  study. 
The  1980  year  group  was  chosen  because  of  the  rather  random  assortment 
of  storms  which  developed  during  that  year  and  the  author's  experience 
and  knowledge  of  them.  Additionally,  the  data  from  Super  Typhoon  Tip 
(1979)  were  included  because  of  the  unique  characteristics  of  this 
typhoon. 

Tropical  cyclone  reconnaissance  missions  are  flown  basically  at  two 
specific  levels.  Low  level  investigative  flights  are  made  during  the 
formative  stages  of  the  system  when  winds  are  below  25  m  s”^ .  These 
flights  are  flown  at  altitudes  of  between  150  and  500  m  above  the  ocean, 
and  sea  level  pressures  are  extrapolated  from  flight  level  D-value  in¬ 
formation.  Once  the  storm  has  definitely  developed  into  a  closed  circu¬ 
lation,  flights  are  conducted  at  the  70  kPa  level  along  radial  legs  both 
into  and  out  of  the  center. 

Data  collection  frequency  also  varies  between  low  level  and  70  kPa 
level  missions.  For  the  low  level  missions,  the  data  of  interest  for 
this  study  were  collected  approximately  every  15  minutes  or  85  km. 
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These  data  included  the  date/time  group,  latitude  and  longitude  of  the 
observation,  flight  level  wind,  sea  level  pressure,  and  surface  wind. 

The  data  collection  frequency  for  70  kPa  missions  was  basically  every 
55  km  within  280  km  of  the  center.  The  data  of  interest  at  this  level 
are  the  date/time  group,  latitude,  longitude,  flight  level  wind,  height 
of  the  70  kPa  surface,  and  the  surface  wind  (if  observable). 

The  initial  flight  path  at  low  level  is  flown  toward  the  forecasted 
center  position  (Fig.  2).  Usually  these  missions  are  flown  before  a 
definite  surface  circulation  has  been  observed.  The  forecasted  position 
usually  is  estimated  from  a  satellite  image  of  the  area  of  disturbed 
weather  and  usually  is  inaccurate  because  of  the  lack  of  cloud  defini¬ 
tion  at  this  stage  of  development.  Even  when  a  circulation  pattern  is 
indicated,  the  circulation  may  be  at  middle  levels  and  may  not  be  evi¬ 
denced  at  the  surface.  When  the  forecasted  position  has  been  reached  by 
the  aircraft,  a  counterclockwise  pattern  is  flown  to  close  off  the 
circulation  (if  any) . 

The  70  kPa  missions  are  flown  in  triangular  crossing  patterns  (Fig. 
3).  Along  the  incoming  radial  leg,  an  observation  is  taken  at  465  km, 
280  km,  and  then  every  55  km  until  the  center  is  fixed.  An  exit  leg  is 
flown  out,  in  the  same  manner,  to  280  km,  and  a  base  leg  is  flown  with 
observations  taken  about  every  130  km.  These  patterns  are  repeated  for 
two  or  three  fixes  of  the  storm's  center.  The  track  between  fixes  re¬ 
quires  about  three  hours  to  complete. 

The  procedures  used  in  obtaining,  sorting  and  processing  the  air¬ 
craft  reconnaissance  data  are  contained  in  the  Appendix.  The  Appendix 
includes  information  on  the  source  of  raw  data,  the  format  in  which  the 


CHAPTER  V 


RESULTS 

In  the  following  discussion,  a  synopsis  will  be  given  of  the  compu¬ 
ter  output  and  its  statistical  significance.  The  applications  of  this 
result  will  be  interpreted  in  terms  of  intensity  and  strength.  Finally, 
the  reliability  of  the  size  measurement  capability  will  be  summarized. 

The  output  of  the  computer  program  for  each  set  of  observations 

with  pressure  data  was  composed  of  ratios  of  the  observed  (average  or 

maximum)  wind  speed  to  the  corresponding  surface  pressure  gradient  or 

70  kPa  height  gradient  between  observations.  Units  of  the  output  were 
2-1-1  2  -1  -1 

in  n  mi  h  mb  or  n  mi  h  f  ,  and  they  were  left  that  way  for 
simplicity  since  the  input  data  were  in  these  terms.  The  values  them¬ 
selves  were  not  as  important  as  their  repeatability,  so  that  for  a  given 
height  or  pressure  gradient  a  realistic  value  of  wind  could  be  assumed 
even  if  winds  were  not  directly  observable  at  the  time  (i.e.,  nighttime 
or  inoperative  wind  sensors).  The  ratios  were  manipulated  statistically 
to  define  mean  relationships  between  wind  and  pressure/height  gradients, 
the  normalcy  of  their  distribution,  and  their  standard  deviation  (or 
variance).  It  was  hoped  that  the  final  results  would  have  a  normal  dis¬ 
tribution,  a  small  variance,  and  that  they  would  support  the  use  of 
empirical  relationships  to  obtain  approximate  wind  speeds  from  observed 
pressure/height  gradients  alone.  The  various  correlations  defined  in 
the  Appendix  were  all  determined  simultaneously  to  ascertain  if  some  of 
the  correlations  were  more  reliable  than  others. 

The  output  values  were  divided  into  groups  which  were  dependent  on 


the  radial  distance  from  the  vortex  center.  Since  data  were  gathered 
approximately  every  55  km  from  the  center,  this  distance  increment  was 
used  as  the  basis  for  the  grouping.  In  doing  so,  the  effect  of  radial 
distance  was  adjusted  in  a  rough  manner  in  the  results.  If  the  midpoint 
between  two  observations  with  pressure  data  lay  within  55  km  of  the 
center,  the  resultant  pressure/wind  relationship  was  placed  in  group  1. 
If  the  midpoint  was  between  56  and  110  km  from  the  center,  the  result 
was  placed  in  group  2,  etc.  The  output  data  also  were  grouped  into  four 
intensity  categories,  using  the  central  pressure  thresholds  that  corre¬ 
late  to  the  transition  point  in  (5)  from  depression  to  tropical  storm, 
tropical  storm  to  typhoon,  and  typhoon  to  super  typhoon.  This  was  done 
to  isolate  those  characteristics  common  to  weaker  storms  from  those 
common  to  stronger  storms. 

A  statistical  analysis  was  completed  to  evaluate  the  distribution 
of  the  output.  Figure  4  compares  the  distributions  of  the  ratios  of 
maximum  surface  wind/sea  level  pressure  gradient  within  55  km  of  the 
center  for  the  four  storm  intensity  categories  (depression,  tropical 
storm,  typhoon,  and  super  typhoon).  The  stem  and  leaf  plots  of  Fig.  4 
approximate  histograms  of  the  output  values  for  the  above  named  distri¬ 
bution  (Koopmans,  1981).  The  box  plots  show  how  the  distribution  com¬ 
pares  to  a  normal  distribution.  Significant  outliers  and  positive  skew¬ 
ness  indicate  that  the  results  are  not  a  normal  distribution.  This  is 
verified  also  by  the  large  differences  between  mean  and  median  values. 

Table  1  compares  the  mean  and  median  values  of  the  distributions 
(similar  to  those  of  Fig.  4)  for  the  ratio  of  average  surface  wind  to 
sea  level  pressure  gradient  and  the  ratio  of  maximum  surface  wind  to  sea 
level  pressure  gradient.  As  an  example,  the  mean  and  median  values  for 
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Fig.  4.  Stem  and  leaf  plots  and  box  plots  for  the  maximum 

surface  wind/sea  level  pressure  gradient  computations 
within  55  km  of  the  center  for  (a)  depressions,  (b) 
tropical  storms,  (c)  typhoons,  and  (d)  super  typhoons. 
The  mild  and  extreme  outliers  of  the  distribution  are 
signified  by  0  and  *,  respectively,  on  the  box  plot. 
Units  in  n  mi 2  h"i  mb-^.  (See  text). 
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Table  1.  Comparison  of  the  mean  (X)  and  median  (M)  values  for 
wind/pressure  relationships  categori2ed  by  intensity  of  and  distance 
from  the  vortex  center.  The  relationships  are  average  surface  wind/sea 
level  pressure  gradient  (ASW/SLPG)  and  maximum  surface  wind/sea  level 
pressure  gradient  (MSW/SLPG).  Units  are  n  mi^  h“^  mb“^. 


Intensity 

Distance  (km) 

ASW/SLPG 

MSW/SLPG 

X 

M 

X 

M 

0-55 

223 

178 

237 

200 

56  -  no 

371 

336 

452 

Depression 

111  -  165 

612 

561 

741 

640 

166  -  220 

436 

493 

502 

221  -  280 

577 

584 

651 

651 

0-55 

177 

153 

190 

165 

56  -  no 

333 

281 

383 

299 

Tropical  Storm 

111  -  165 

504 

419 

548 

452 

166  -  220 

470 

404 

536 

456 

221  -  280 

653 

546 

743 

0  -  55 

94 

83 

100 

91 

56  -  no 

238 

185 

254 

198 

Typhoon 

111  -  165 

382 

314 

418 

341 

166  -  220 

518 

425 

555 

448 

221  -  280 

500 

416 

311 

416 

0  -  55 

37 

34 

39 

35 

56  -  no 

218 

195 

242 

198 

Super  Typhoon 

111  -  165 

283 

256 

293 

275 

166  -  220 

534 

411 

558 

448 

221  -  280 

706 

567 

546 

628 

the  0-55  km  distance  range  in  Table  1  correspond  to  the  mean  and  median 
values  of  the  stem  and  leaf  plot  in  Fig.  4a.  Thus,  the  mean  ratio  of 
the  average  surface  wind  (ASW)  to  the  sea  level  pressure  gradient  (SLPG) 
is  223,  whereas  the  median  ratio  value  is  178.  The  other  correlations 
that  were  tabulated  by  the  computer  (i.e.,  70  kPa  wind/70  kPa  height 
gradient  and  surface  wind/70  kPa  height  gradient,  etc.)  exhibited  simi¬ 
lar  tendencies  and  are  not  shown.  The  effects  of  the  mild  or  extreme 
outliers  skew  the  means  considerably,  widen  the  difference  between  mean 
and  median,  destroy  distribution  normalcy,  and  limit  the  data's  useful¬ 
ness  in  determining  empirical  relationships. 

Figure  4a  indicates  the  presence  of  significant  negative  outliers. 
These  are  generally  a  result  of  inaccurate  pressure  gradient  computa¬ 
tions,  incorrect  wind  measurements  (or  estimates),  or  strong  local  wind 
effects.  Pressure  gradient  inaccuracies  can  be  attributed  to  the 
effects  of  rounding  error,  poor  navigational  positioning,  pressure 
measurement  error,  or  possibly  even  an  invalid  sea  level  pressure  extra¬ 
polation  assumption.  Local  wind  effects,  such  as  a  thunderstorm  out¬ 
flow,  also  may  be  contrary  to  the  assumed  wind  of  the  observed  pressure/ 
height  gradient.  The  negative  relationship  commonly  occurs  in  weak 
pressure  gradient  situations  (0.1  -  0.2  kPa  differences  between  observa¬ 
tions),  in  broad  light  and  variable  wind  areas,  or  where  the  flight  path 
is  nearly  parallel  to  the  height  contours  or  isobars  of  the  pressure 
field.  The  same  inaccuracies  also  may  account  for  some  of  the  positive 
outliers.  More  restrictive  requirements  covering  the  normalcy  of  the 
flight  path  to  the  pressure  field  may  be  necessary  to  prevent  these  data 
from  affecting  the  results.  The  wind  measurement  itself  is  another 


factor  in  the  variability  of  the  outcome.  The  ability  to  determine  cor¬ 
rectly  the  average  or  maximum  wind  speeds  is  subject  to  question  when 
only  two  or  three  wind  measurements  are  taken  along  a  55  km  track.  Add¬ 
ed  to  this  are  the  factors  of  Doppler  radar  attenuation  in  heavy  precip¬ 
itation  and  unobservable  surface  winds  during  undercast  cloud  condi¬ 
tions  . 


A.  Intensity 


The  gradient  winds  were  computed  using  the  equation. 


(6) 


where  is  the  gradient  wind,  f  is  the  Coriolis  term,  r  is  the  radial 
distance  from  the  vortex  center  to  the  midpoint  between  observations, 
the  specific  gas  constant,  T^  is  the  virtual  temperature  (assumed  to  be 
a  constant  29°C),  dp  the  change  in  pressure  between  observations,  p  is 
the  assumed  pressure  at  the  midpoint,  and  dr  the  distance  between  obser¬ 
vations  along  the  same  radial  line.  For  small  values  of  r  (less  than 
55  km)  at  low  latitudes,  the  term  is  negligible  such  that 


where  is  the  cyclostrophic  wind. 

Tables  2-5  contain  comparison  summaries  of  randomly  selected  storms 


in  all  four  intensity  categories.  The  tables  compare  the  computations 
of  surface  winds  from  the  computer  program  (V -  ),  the  maximum  wind 
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Table  2.  Comparison  of  the  computed  and  observed  maximum  surface 
winds  within  55  km  of  the  center  of  a  tropical  depression.  Units  are  in 


m  s 


-1 
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Winds  as 
V.u. 


Table  3.  As  Table  2,  except  for  tropical  storms, 


as  computed  by  different  methods 


different  methods 


Observed 
maximum 
surface  wind 


23 

33 

20 

26 

19 

25 

16 

9 

16 

23 

16 

15 

20 

24 

20 

24 

21 

29 

22 

22 

Observed 
maximum 
surface  wind 


Winds  as 

computed  by  different 

means 

Observed 

'^cbs 

Va. 

maximum 
surface  wind 

•  •  •  •  , 

..  59 

71 

59 

63 

72 

71 

60 

66 

61 

75 

63 

67 

74 

75 

64 

66 

85 

83 

66 

69 

61 

83 

67 

63 

U^-il 

114 

82 

61 

74 

77 

82 

63 

67 

55 

79 

64 

69 

101 

73 

61 

62 

71 

73 

67 

64 

80 

75 

69 

68 

• 

73 

75 

68 

71 

based  on  the  Atkinson  and  Holliday  equation  (V^^),  the  gradient  wind 
(Vgj,),  and  the  observed  maximum  surface  wind.  We  must  assume  that  the 
maximum  winds  will  be  found  within  55  km  of  the  center,  although  this 
may  not  be  the  case  for  depressions  and  other  disorganized  storms.  In 
comparing  the  three  methods,  the  variability  of  using  the  data  empiri¬ 
cally  is  evident,  especially  for  typhoons  and  super  typhoons 

(Tables  4  and  5),  and  it  is  more  realistic  to  assume  that  the  gradient 
wind  and  the  Atkinson  and  Holliday  methods  are  better. 

Figure  5  is  a  plot  of  the  average  differences  between  and 
for  the  four  intensity  categories.  As  the  storm  intensifies,  in¬ 
creasingly  underestimates  the  maximum  wind  (if  we  assume  to  be 
correct).  This  can  be  explained  by  the  classic  pressure  configuration 
within  the  eyes  of  tropical  cyclones.  Even  without  an  organized  eye 
wall,  the  pressure  profile  cannot  be  approximated  by  the  gradient  wind 
relationship  once  inside  the  maximum  wind  band.  Many  theoretical  wind 
approximations  of  the  tropical  cyclone  have  assumed  solid  rotation  with¬ 
in  the  radius  of  maximum  wind.  In  other  words,  the  winds  are  expected 
to  decrease  to  zero  in  a  linear  fashion  from  the  point  of  maximum  wind 
to  an  infinitely  small  calm  center.  In  reality,  this  is  not  the  case. 
Most  storms  show  a  broad  minimum  pressure  center  with  light  and  variable 
or  calm  winds.  These  minimum  pressure  centers  commonly  have  radii  of 
10-20  km.  In  Fig.  6,  the  actual  minimum  pressure  center  of  radius  rp  is 
characterized  by  virtually  no  pressure  gradient  or  measurable  wind.  It 
is,  therefore,  useless  to  include  this  distance  in  the  innermost  estima¬ 
tion  of  Vg^.  If,  however,  rp  were  used  to  adjust  the  distance  over 
which  the  pressure  gradient  is  measured  (change  r^  to  r^)  and  to  adjust 
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Average  differences  between  V  .  and  V  for  the  four  tropical 
cyclone  intensity  categories.®  ^ 


Fig,  6.  Adjustment  of  the  distance  over  which  pressure  gradient  is 

measured  (r^  to  r  ')  and  the  radius  of  measurement  (r  to  r  ). 
The  minimum  sea  level  pressure  radius  (r  )  represents  the  ^ 
point  on  the  pressure  curve  after  which  a  quasi-1 inear 
approximation  can  be  made. 


the  radial  distance  to  the  midpoint  for  which  this  gradient  applies 
(change  r  to  r,),  then  by  substituting  r,  for  r  and  r'  for  dr  in  (7), 
the  differences  between  and  would  be  approximately  offset.  As 
an  example,  a  super  typhoon  at  15°  latitude  with  minimum  sea  level 
pressure  of  89  kPa  and  a  radial  distance  of  55  km  (r^)  to  an  observatipn 
with  a  pressure  rise  (from  the  central  pressure)  of  8.6  kPa  (dp  =  8.6 
kPa)  would  require  a  minimum  pressure  center  radius  (rp)  of  9  km  to  off¬ 
set  the  12  m  s~^  difference  between  and  the  original  This  is  a 

realistic  value  for  storms  of  this  intensity. 

Figure  7  shows  that  not  only  is  ^(Vg^  -  Vg^)  /  N  >  0  for  typhoons 
and  super  typhoons,  but  the  same  relation  holds  true  for  the  average 
differences  between  and  the  maximum  surface  winds  observed  (estima¬ 
ted)  by  aircraft  observers  (V^).  Although  ^(Vg^  -  Vg^)  /  N  and 
r(Vah  ■  Vq)  /  N  exhibit  similar  trends,  the  degree  to  which  is  less 
than  Vg^  is  accounted  for  in  a  different  way.  During  the  depression 
stage,  there  is  generally  less  cloudiness  below  observational  altitude 
that  will  interfere  with  the  observation  of  the  maximum  surface  wind. 
Additionally,  flights  are  frequently  flown  at  low  levels  during  the 
depression  stage,  which  allows  for  more  accurate  wind  estimations.  As 
storm  strength  increases  and  an  eye  wall  forms,  the  chances  increase 
that  the  maximum  winds  will  be  missed  while  the  aircraft  is  in  heavy 
precipitation  or  the  eye  wall  cloud.  Estimations  of  surface  winds  above 
67  m  s”^  are  further  complicated  by  the  inability  to  distinguish  changes 
when  the  surface  is  obscured  by  spray.  With  the  adjustment  for  system 
motion,  five  of  the  thirteen  observed  winds  in  Table  5  are  above  67  m 
s"^ ;  however,  only  two  of  the  original  wind  estimates  (V  )  were  for 


I  DEPREJ/^  I  TROPICAL  j  TYPHOON  |  SUPER  TYPHOON 

SION*  STORM 

STORM  INTENSITY 


As  Fig.  5,  except  between  V  .  and  V  ,  V  .  and  V.,  as  well  as 
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winds  in  excess  of  67  m  s'^ . 

Figure  7  also  shows  a  similar  comparison  for  instrumental ly  observ¬ 
ed  maximum  winds  at  flight  level  (V^).  The  same  tendency  occurs  with 
increasing  storm  intensity  as  was  the  case  for  the  observations  of  sur¬ 
face  wind.  This  is  due  to  the  general  increase  in  heavy  precipitation 
with  the  more  intense  storms.  Additionally,  Doppler  radar  systems  often 
reach  their  maximum  drift  angle  measurement  capability  during  penetra¬ 
tion  of  the  eye  wall,  preventing  instrumental  measurement  of  the  maximum 
winds. 

We  must  either  assume  that  overestimates  the  maximum  wind  or 
consider  both  measured  and  estimated  winds  from  aircraft  as  unreliable. 
The  latter  hypothesis  is  preferred  in  light  of  the  high  variability  of 
the  computed  wind  data  of  this  study  and  the  known  problems  of  accurate¬ 
ly  observing  winds  from  aircraft.  However,  with  an  accurate  approxima¬ 
tion  of  the  radius  of  minimum  pressure  change  (r^)  within  tropical 
cyclone  centers  and  accurate  position  reports,  a  meaningful  intensity 
approximation  is  attainable  from  the  aircraft  pressure  data. 

B.  Strength 

Figures  8-11  are  wind  profiles  of  selected  radial  legs  flown  into 
or  out  of  tropical  cyclones  of  all  four  intensity  categories.  All  winds 
in  the  figures  have  had  the  motion  of  the  system  subtracted.  As  would 
be  expected  for  tropical  depressions  (Fig.  8),  there  is  no  classic  wind 
or  pressure  relationship  such  as  is  common  for  much  stronger  storms. 

The  profile  for  Tropical  Storm  Forest  is  for  the  period  approximately  18 
hours  before  he  was  classified  as  a  tropical  storm.  The  closest 
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observation  to  the  center  indicated  that  there  were  surface  winds  of  15 
m  s”^  at  96  km  from  the  center,  but  from  that  point  on  into  the  center 
there  is  no  information.  The  pressure  gradient  information  seems  to 
support  the  observed  wind,  and  a  secondary  maximum  of  at  174  km 
(21  m  s'^)  would  seem  to  justify  the  observed  maximum  surface  winds  of 
20  m  s'^  at  152  km.  The  minimum  sea  level  pressure  of  99.9  kPa  would 
support  a  maximum  wind  of  only  16  m  s”^  using  (5),  and  without  knowing 
of  the  maximum  outside  100  km  there  would  be  no  reason  to  suspect 
the  existence  of  stronger  winds  at  that  range.  But  successive  tracks 
through  the  area  observed  the  strongest  winds  between  100  and  130  km 
from  the  center. 

The  profile  of  what  would  become  Typhoon  Tip  is  also  included  in 

Fig.  8.  The  gradient  wind  would  seem  to  indicate  that  winds  of  10  m  s”^ 

were  justified  outside  of  80  km  and  would  also  seem  to  indicate  that 
winds  would  pick  back  up  again  within  70  km  of  the  storm.  This  was  not 
to  be  the  case.  To  help  explain  the  reason  for  this,  we  must  bring. in 
the  concept  of  vertical  tilt  of  the  storm.  In  order  for  to  be  com¬ 
puted  in  (6),  the  surface  pressures  were  used  to  determine  dp.  However, 
the  surface  pressures  themselves  were  extrapolated  from  70  kPa  heights, 
which  are  related  to  the  70  kPa  center  but  not  directly  related  to  the 
surface  center.  Only  when  surface  and  70  kPa  centers  are  vertically 
aligned  can  we  assume  that  extrapolated  sea  level  pressure  from  70  kPa 

data  will  duplicate  actual  sea  level  pressures  at  the  surface.  If  the 

surface  center  is  displaced  in  the  horizontal  from  the  70  kPa  center, 
the  extrapolated  sea  level  pressures  will  be  unlikely  to  approximate  the 
actual  surface  values  and  estimating  dp  becomes  unreliable  even  though 
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dr  between  the  observation  point  and  both  centers  is  accurately  deter¬ 
mined.  In  this  particular  case,  the  vertical  displacement  of  the  storm 
was  judged  to  be  70-80  km  between  sea  level  and  70  kPa  (Huntley  and 
Diercks,  1981).  This  complicated  the  accuracy  of  sea  level  pressure 
extrapolations  and  subsequent  computations  of  Vg^  close  to  the  center. 

The  profile  for  Tropical  Storm  Kim  is  shown  for  successive  tracks 
in  two  separate  sectors  (Fig.  9).  Following  the  route  of  flight  on  the 
inbound  track,  it  is  interesting  to  note  that  a  relatively  strong  wind 
zone  exists  at  around  200  km.  The  value  of  Vg^  computed  at  235  km  be¬ 
tween  the  observations  at  198  and  272  km  indicates  that  winds  should  in¬ 
crease,  although  the  observed  winds  do  not  increase  as  much  as  predic¬ 
ted.  The  Vg^  computed  for  181  km  is  significantly  higher  than  the  ob¬ 
served  wind,  which  actually  decreases  by  164  km  as  compared  with  the 
value  at  198  km.  Considering  that  the  navigational  positioning  error 
may  be  11  km  or  more  (positions  are  rounded  to  the  nearest  tenth  of  a 
degree  in  latitude  and  longitude),  an  adjustment  to  the  observational 
position  from  164  km  to  153  km  would  decrease  Vg^  at  181  km  by  4  m  s  ^ . 

This  value  of  V  would  still  appear  to  overestimate  the  observed  wind 
gr 

but  would  indicate  that  a  slight  decrease  in  wind  speed  (as  observed) 
was  likely.  The  increase  in  Vg^  at  72  km  is  reflected  in  a  slight  in¬ 
crease  in  surface  wind  at  48  km  without  a  corresponding  increase  in 
flight  level  winds.  This  may  be  due  to  vertical  tilt  effects  since  the 
inbound  observations  were  taken  from  a  position  4-5  km  closer  to  the 
70  kPa  center  than  to  the  surface  center.  On  the  outbound  leg,  the  ob¬ 
served  maximum  winds  that  occurred  at  174  km  appear  justified  by  Vg^  if 
we  assume  that  the  pressure  gradient  increase  takes  place  right  at  the 
174  km  observation.  If  we  adjust  the  174  km  position  by  11  km  to  163 
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km,  the  computed  at  139  km  goes  to  16  m  s  and  at  200  km  falls 
to  34  m  s"^ .  This  would  appear  more  reasonable,  but  the  tendency  of 
to  overestimate  observed  wind  values  remains. 

The  more  classic  profile  appears  for  typhoons  (Fig.  10).  The  in¬ 
bound  track  on  Typhoon  Lex  shows  that  does  a  fair  job  of  approximat¬ 
ing  flight  level  winds  but  continues  to  exceed  the  estimates  of  surface 
wind.  The  inbound  plot  on  Typhoon  Wynne  portrays  a  very  strong  storm. 
Typhoon  force  winds  appear  to  extend  to  nearly  200  km.  As  this  storm 
moves  closer  to  a  strong  subtropical  ridge  over  the  Asian  continent  and 
North  Pacific,  a  strengthened  gradient  extends  outward  to  the  north  of 
the  storm,  and  the  storm  seems  to  strengthen  even  though  it  is  now  past 
the  super  typhoon  stage.  The  computation  of  above  50  m  s"^  at  256 
km  does  not  appear  to  be  justified  and  may  be  associated  with  an  obser¬ 
vational  error. 

Figure  11  is  the  plot  of  wind  speeds  along  consecutive  inbound  and 
outbound  tracks  for  Super  Typhoon  Tip.  The  observed  winds  on  the  in¬ 
bound  track  appear  to  support  that  is  too  high.  For  example,  the 
values  of  at  76,  128,  180,  231,  and  354  km  are  all  above  the  traces 
which  connect  the  observed  surface  winds  (V^)  and  instrumental 1y  sensed 
flight  level  winds  (V.j).  In  fact,  this  would  tend  to  indicate  that  for 
the  five  computations  of  outside  of  70  km,  was  greater  than 

by  around  7  m  s”^  and  V.j  by  over  6.5  m  s”\  However,  the  outbound  track 
shows  only  slight  differences  between  and  V^.  Even  though 

maximum  observed  and  instrumental 1y  sensed  winds  are  used  in  the 
figures,  these  winds  may  be  underestimates  of  the  true  maximum  wind 
values  between  observations.  It  was  discussed  that  there  are 


difficulties  in  observing  or  sensing  maximum  winds  with  the  limited  op¬ 
portunities  afforded  by  extensive  cloud  cover  (which  limits  the  view  of 
the  sea  surface)  and  heavy  rainfall  (which  attenuates  Doppler  speed 
measurements).  Super  Typhoon  Tip  had  both  of  these  traits,  and  it  is 
quite  likely  that  one  quadrant  of  the  storm  (inbound  track)  was  not  con¬ 
ducive  to  measuring  maximum  values  while  another  quadrant  was  conducive 
(outbound  track).  Given  that  there  is  a  good  correlation  between 
V^,  and  V^.  on  the  outbound  track  and  since  certainly  does  not  under¬ 
estimate  maximum  values  of  V„  and  V.,  it  can  be  assumed  that  is  a 

0  1’  gr 

good  indicator  when  used  to  determine  storm  strength. 

A  comparison  of  the  "strength"  of  Super  Typhoons  Tip  (with  p^.  = 

90  kPa)  and  Wynne  (with  p^  =  89  kPa)  is  made  in  Fig.  12  for  both  and 
Vgr-  The  comparison  is  made  while  the  central  pressure  values  (p^.)  of 
both  storms  are  within  1  kPa  of  each  other,  not  necessarily  when  they 
are  the  lowest  value.  Both  V„  and  indicate  that,  at  this  point. 
Super  Typhoon  Wynne  is  more  "intense"  than  Tip  at  the  maximum  wind  band 
(as  p^  would  indicate).  However,  Tip  exhibits  greater  "strength"  out¬ 
side  of  110  km.  As  Typhoon  Wynne  was  downgraded  from  super  typhoon  to 
typhoon  (with  p^  =  92.5  kPa),  she  continued  to  strengthen  outside  of 
100  km  (Fig.  12).  Even  though  the  value  of  at  256  km  does  not  seem 
justified,  as  discussed  earlier,  the  environmental  impact  is  evident 
when  comparing  the  strength  of  these  two  storms,  and  central  pressure 
seems  to  have  little  impact. 

If  not  used  quantitatively,  the  qualitative  use  of  Vg^  can  approxi¬ 
mate  the  characteristic  strength  of  a  tropical  cyclone.  Also,  it  may  be 
used  to  determine  the  wind  profiles  of  disorganized  systems  with 


unusually  high  winds  occurring  well  separated  from  the  center.  In  order 
to  do  this  accurately,  a  finer  mesh  of  pressure  observations  will  need 
to  be  constructed  with  accurate  instrumentation  and  precise  navigation. 
This  is  not  feasible  with  present  operational  methods  of  typhoon 
reconnaissance. 

C.  Size 

The  storm  size  was  evaluated  for  all  tropical  cyclones  in  the  data 
set.  Tables  6  and  7  contain  a  comparison  of  the  various  techniques  used 
in  estimating  the  size  of  four  storms.  The  computer  program's  computa¬ 
tion  of  size  is  compared  to  various  manual  techniques  using  the  synoptic 
scale  analysis  that  is  available  over  the  ocean.  The  primary  synoptic 
means  of  determining  size  is  the  AROCI  technique  discussed  in  Chapter 
III.  The  Daily  Weather  Maps  of  the  Japan  Meteorological  Agency  were 
used  for  the  synoptic  scale  surface  analysis.  On  this  scale,  which  is 
the  only  scale  available  over  data  sparse  oceans,  the  outer  closed  iso¬ 
bar  was  nearly  always  analyzed  as  100.8  kPa.  Therefore,  an  accurate 
estimate  of  the  actual  outermost  closed  isobar  value  was  not  always 
possible.  Yet,  using  100.8  kPa,  the  AROCI  values  are  listed  in  the 
third  column  of  the  tables. 

Often,  the  100.8  kPa  isobar  was  extremely  elongated.  This  would 
tend  to  yield  unrealistically  large  values  for  the  AROCI.  When  this  was 
the  case,  the  next  lower  symmetric  isobar  at  pressure  p  was  subjective¬ 
ly  approximated  by  interpolation  between  analyzed  values  to  compute  the 
average  radius  of  the  outermost  closed  symmetrical  isobar  (AROCSI). 

These  values  were  preferred  over  the  AROCI  and  also  are  listed  in  the 
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Table  6.  Comparison  table  of  storm  size  computation  methods  for 
tropical  depressions  and  tropical  storms.  The  table  compares  the  pro¬ 
gram's  computed  radius  (r'),  the  AROCI  using  100.8  kPa,  the  AROCI  using 
a  variable  pressure  (p'h  and  the  average  radius  of  the  outer  closed 
symmetrical  isobar  (AROCSI)  at  pressure  p. 


Row 

no. 

(km) 

AROCI  (100.8  kPa)  AROCSI/p 

(km)  (km/kPa) 

AROCI/p' 

(km/kPa) 

NAME 

1. 

306 

363 

DEPRESSION 

165/100.5 

DOM 

2. 

461 

265 

- 

552/101.1 

DOM 

3. 

326 

219 

- 

363/100.9 

FORREST 

4. 

189 

222 

- 

400/101.0 

FORREST 

5. 

265 

222 

- 

- 

FORREST 

6. 

294 

267 

- 

457/101.0 

FORREST 

7. 

889 

544 

- 

193/100.6 

TDM 

1. 

235 

372 

TROPICAL  STORM 

268/100.6 

DOM 

2. 

237 

417 

- 

354/100.7 

DOM 

3. 

415 

698 

524/100.6 

- 

DOM 

4. 

345 

578 

- 

- 

DOM 

5. 

237 

- 

- 

469/100.7 

DOM 

6. 

402 

- 

- 

485/100.7 

DOM 

7. 

180 

244 

- 

• 

FORREST 

8. 

463 

270 

• 

• 

FORREST 

9. 

222 

217 

- 

531/101.0 

FORREST 

10, 

146 

137 

* 

- 

FORREST 

11. 

293 

407 

• 

- 

TIP 

12, 

443 

1141 

420/100.4 

913/100.7 

TIP 

13. 

396 

1151 

354/100.4 

1151+/100.9 

TIP 

14. 

243 

930 

257/100.4 

330/100.5 

TIP 

15. 

754 

931 

293/100.4 

- 

TIP 

•  .1 
%  •'  4 
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tables  when  extreme  elongation  of  the  100.8  kPa  isobar  occurred. 

As  mentioned  earlier,  the  program  automatically  assumes  that  the 
outer  limit  of  the  storm  has  been  reached  when  flight  level  and  surface 
winds  drop  below  13  and  11  m  s'^  respectively.  Initially,  the  pressure 
at  the  outer  limit  of  the  storm  (p')  is  assumed  to  be  100.8  kPa;  how¬ 
ever,  when  the  above  wind  limits  are  reached  at  lower  pressures,  the 

program  uses  this  new  value  of  p'  in  all  computations.  Conversely,  when 
winds  at  higher  pressure  values  (p  >  p')  are  above  18  and  13  m  s  ^  re¬ 
spectively,  p'  is  changed  upward.  The  synoptic  analysis  of  the  AROCI 
using  p'  is  then  estimated  and  given  preference  over  the  synoptic  analy¬ 
sis  of  the  AROCI  using  100.8  kPa  and  that  of  the  AROCSI  using  p  as  the 
best  synoptic  means  of  determining  size. 

In  the  tables,  the  preferred  synoptic  means  of  analyzing  the  AROCI 
is  listed  farther  to  the  right.  In  the  case  of  depressions,  the  dif¬ 
ferences  between  the  computed  radius  (r')  of  p'  and  the  analyzed  AROCI 

values  for  p'  suggest  that  the  precision  of  the  AROCI  technique  on  the 
synoptic  scale  is  questionable.  It  is,  therefore,  not  advisable  to  re¬ 
gard  the  preferred  synoptic  means  of  computing  the  AROCI  as  absolute  in 
defining  the  size  of  the  storm. 

Generally  speaking,  the  r'  values  that  are  given  in  the  tables  com¬ 
pare  favorably  with  the  preferred  AROCI  values  when  p'  is  measured 
directly  on  one  of  the  inbound  or  outbound  tracks  or  when  r'  is  inter¬ 
polated  using  (A-7).  However,  in  those  cases  where  data  were  not  taken 
at  sufficient  distance  from  the  storm  for  interpolation,  r'  tends  to 
blow  up  to  unacceptable  values.  This  can  be  seen  for  the  last  depres¬ 
sion  (row  no.  7)  and  tropical  storm  (row  no.  15)  listed  in  Table  6  and 


for  several  typhoon  and  super  typhoon  cases  in  Table  7. 

A  search  was  made  for  ship  reports  to  supplement  the  aircraft  data. 
Reconnaissance  data  for  missions  flown  near  synoptic  report  times  on 
TD  14  and  Typhoon  Tip  were  supplemented  with  nearby  synoptic  ship  re¬ 
ports.  The  ship  report  position  and  pressure  were  added  to  the  data  so 
that  r'  could  be  interpolated  instead  of  extrapolated.  The  corrected 
values  of  r'  were  computed  with  (A-7),  using  the  computer's  output  of  A 
and  B  scaling  parameters,  and  are  listed  in  Table  8  along  with  the 
corresponding  improvement.  Improvement  was  calculated  by  determining 
the  absolute  difference  between  r'  and  the  preferred  AROCI  technique  and 
subtracting  the  absolute  difference  between  the  corrected  r'  and  the 
preferred  AROCI  technique.  This  is  only  meaningful  if  we  assume  that 
the  preferred  AROCI  is  the  best  standard  for  comparison.  General  im¬ 
provement  was  noted  for  most  cases,  but  some  cases  indicated  degrada¬ 
tion.  However,  examination  of  the  initial  reconnaissance  data  for  the 
first  typhoon  listed  in  Table  8  reveals  that  a  100.8  kPa  surface  pres¬ 
sure  was  extrapolated  from  the  70  kPa  level  at  365  km  from  the  center 
(versus  a  corrected  r'  value  of  369  km  interpolated  by  the  program  with 
the  ship  report  included).  This  information  was  not  used  in  the  origi¬ 
nal  computation  of  r'  because  the  scaling  parameters  had  not  yet  been 
computed.  However,  it  does  show  that  the  AROCI  technique,  in  itself,  is 
subject  to  considerable  error  without  the  use  of  aircraft  and  ship  in¬ 
formation.  Very  little  ship  pressure  information  was  incorporated  into 
the  Japanese  analysis  even  though  ship  wind  reports  were  used. 

When  the  contents  of  Table  8  are  plotted  for  Typhoon  Tip  in  its 
typhoon  and  super  typhoon  stages  (Fig.  13),  the  extreme  variation  of  the 
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Table  8.  Improvement  of  r'  computation  with  the  addition  of  synop¬ 
tic  ship  reports. 


r' 

(km) 


Preferred 

AROCI/p 

(km/kPa) 


Corrected  r' 
(km) 


Improvement 

(km) 


Name  (row  no.) 


DEPRESSION 


889 

193/100.6 

269 

TROPICAL  STORM 

620 

TD14 

(7) 

415 

524/100.6 

569 

64 

DOM 

(3) 

345 

578/100.8 

409 

TYPHOON 

64 

DOM 

(4) 

843 

861/100.8 

369 

-474 

TIP 

(8) 

454 

944/100.8 

746 

292 

TIP 

(9) 

1044/100.8 

1135 

-  31 

TIP 

(10) 

2735 

1030/100.8 

748 

1423 

TIP 

(11) 

959 

754/100.4 

725 

176 

TIP 

(12) 

928 

839/100.6 

1050 

-122 

TIP 

(13) 

985/100.8 

815 

921 

TIP 

(14) 

846/100.8 

468 

-213 

TIP 

(15) 

972 

937/100.8 

502 

SUPER  TYPHOON 

-400 

TIP 

(16) 

1870 

1050/100.8 

730 

500 

TIP 

(1) 

639 

1161/100.8 

1226 

457 

TIP 

(2) 

572 

1143/100.8 

1065 

493 

TIP 

(3) 

1470 

1222/100.8 

928 

-  46 

TIP 

(4) 

725 

1176/100.8 

1113 

388 

TIP 

(5) 

Note:  Numbers  in  parentheses  in  the  right  column  refer  to  row 
numbers  in  Table  6  (for  depressions  and  tropical  storms) 
and  Table  7  (for  tvohoons  and  super  typhoons). 


RADIUS  (km) 
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TIME  (days) 


Fig.  13.  Comparison  of  the  values  of  r'  (solid  circles),  corrected  r' 

(solid  triangles),  and  preferred  AROCI  technique  (open  circles) 
for  those  cases  in  Super  Typhoon  Tip's  life  cycle  for  which  p’ 
and  r'  were  extrapolated. 


uncorrected  computer  values  of  r'  is  evident.  With  the  addition  of  only 
one  ship  report  per  mission,  the  size  increase  of  the  storm  is  very 
evident  and  a  considerable  smoothing  is  realized.  With  the  addition  of 
more  supplemental  surface  reports  into  the  model  (either  aircraft,  ship, 
or  surface  observation)  the  model  will  be  able  to  portray  storm  size  in 
real  time. 
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CHAPTER  VI 


CONCLUSIONS 

The  ability  of  routine  aircraft  reconnaissance  data  to  depict 
accurately  the  intensity,  strength,  and  size  characteristics  of  western 
North  Pacific  tropical  cyclones  was  studied  for  representative  tropical 
cyclones  of  1979-1980.  Based  upon  a  statistical  analysis  of  the  useful¬ 
ness  of  aircraft  wind  observations  and  a  comparison  with  computed  gradi¬ 
ent  winds  (using  the  pressure  data),  the  following  conclusions  are  made. 

When  observed  winds  (instrumentally  measured  and  visually  estima¬ 
ted)  are  correlated  with  the  pressure  data,  the  results  show  a  high  de¬ 
gree  of  variance.  This  is  attributed  to  various  instrumental,  observa¬ 
tional,  and  navigational  shortcomings.  These  shortcomings  prevent  the 
data  from  being  used  to  develop  empirical  relationships,  since  the  re¬ 
sultant  data  distributions  do  not  resemble  normal  distributions.  Cer¬ 
tain  steps  may  be  taken  to  "sanitize"  the  data  of  some  of  the  observa¬ 
tional  shortcomings.  For  example,  data  could  be  sorted  or  classified  to 
include  those  observations  which  are  taken  along  tracks  roughly  normal 
to  the  pressure  gradient,  but  to  exclude  observations  along  tracks 
approximately  parallel  to  the  pressure  gradient.  The  exclusion  of  these 
latter  data  would  more  than  likely  negate  some  of  the  problems  which 
account  for  the  extreme  outliers  in  the  wind/pressure  gradient  distribu¬ 
tions.  However,  this  would  still  leave  the  problems  of  inaccurate  wind 
measurement/estimation  and  navigation.  With  this  in  mind,  it  is  assumed 
that  the  means  of  collecting  pressure  data  are  far  more  reliable  than 
quantitatively  determining  wind  speed;  this  assumption  favors  the 


computation  of  gradient  winds  as  an  estimate  of  the  local  wind. 

The  locally  computed  gradient  wind  may  be  used  to  approximate  the 
intensity  of  tropical  cyclones.  For  classically  developed  typhoons, 
which  exhibit  a  sharp  pressure  drop  and  maximum  winds  within  55  km  of 
the  center,  the  gradient  (or  cyclostrophic)  wind,  that  is  computed  using 
the  pressure  gradient  over  this  distance,  underestimates  the  theoretical 
maximum  winds  determined  by  (5).  However,  when  the  radius  of  the 
storm's  minimum  pressure  center  is  known,  the  corresponding  adjustments 
to  pressure  gradient  and  the  radius  of  curvature  terms  in  the  gradient 
(and  cyclostrophic)  wind  equation  lessens  the  differences  between 
methods.  With  improved  ability  to  determine  the  radius  of  negligible 
pressure  change  within  the  eye  and  with  shorter  distances  over  which  to 
measure  pressure  gradients,  the  gradient  wind  method  may  be  a  better 
method  for  determining  intensity  characteristics  of  different  storms 
with  approximately  equal  central  pressures.  The  main  shortcoming  con¬ 
tinues  to  be  navigational  error.  Gradient  wind  errors  of  as  much  as 
6  m  s"^  are  common  for  the  typical  round  off  of  positions  to  the  nearest 
tenth  of  a  degree.  This  must  be  improved  upon  to  be  useful. 

The  disorganized  storms  which  do  not  exhibit  classical  pressure  and 
wind  profiles  also  may  be  evaluated  for  intensity  using  the  local  gradi¬ 
ent  wind.  A  review  of  certain  cases  shows  that  some  disorganized  storms 
exhibit  maximum  wind  bands  at  considerable  distances  from  their  centers. 
These  trends  are  reflected  weakly  by  the  gradient  wind  profiles.  With 
higher  navigational  precision  and  a  denser  data  collection  frequency, 
these  features  may  be  better  defined. 

The  computation  of  gradient  winds  at  approximate  55  km  intervals 


also  may  be  used  to  estimate  general  strength  characteristics.  With  the 
current  observation  interval  as  large  as  55  km  and  with  the  current 
positioning  error  as  large  as  9-18  km,  exact  profiles  are  not  obtain¬ 
able;  however,  general  strength  patterns  may  be  subjectively  made. 

Actual  wind  measurements/estimates  tend  to  be  less  than  the  computed 
gradient  winds;  however,  the  actual  maximum  winds  between  observations 
are  likely  to  be  missed  due  to  a  low  sampling  rate  or  measurement  in¬ 
accuracy. 

The  fitting  of  known  pressure  data  to  a  theoretical  tropical 
cyclone  pressure  profile  is  used  to  approximate  storm  size.  The  short¬ 
comings  of  trying  to  fit  data  from  nonclassical  or  disorganized  storms 
(i.e.,  depressions,  tropical  storms,  etc.)  to  a  classical  profile  are 
at  least  partially  offset  by  internal  program  wind  thresholds,  which  act 
as  a  check  and  balance  system.  However,  large  variations  in  computed 
size  values  point  out  the  need  for  the  following: 

1)  Since  extrapolation  of  radial  size  results  in  the  most  varia¬ 
tion,  data  should  be  collected  at  the  70  kPa  level  at  least  as  far  out 
from  the  center  as  the  storm's  radial  size  (r').  This  would  avoid  the 
need  for  extrapolation  of  the  radial  size  using  (A-7). 

2)  If  extrapolation  is  necessary,  subjective  determination  of  p' 
and  Pj^  may  be  necessary  in  order  to  avoid  excessive  error  in  the  size 
computation.  Small  differences  (on  the  order  of  0.1  kPa)  in  p'  or  p^ 
may  cause  significant  differences  in  computed  size  under  these  circum¬ 
stances. 

3)  Supplementing  aircraft  data  with  as  little  as  one  ship  report 
or  surface  report  provides  an  excellent  means  of  avoiding  extrapolation. 


The  innermost  ship  or  surface  observation  outside  of  r'  with  winds  less 
than  13  m  s"^  would  be  a  good  addition  to  the  aircraft  data  base,  pro¬ 
vided  the  times  of  the  observations  are  approximately  the  same.  The 
ideal  situation  would  be  that  the  computations  in  all  four  quadrants  of 
the  storm  be  averaged  similarly  to  the  AROCI  technique. 

The  operational  aircraft  horizontal  observation  system  has  remained 
essentially  unchanged  for  over  15  years.  These  manual  techniques  of 
data  gathering  do  not  meet  the  needs  of  this  study,  much  less  the  other 
prediction  techniques  that  are  currently  operational.  A  high  density 
horizontal  weather  observation  system  similar  to  the  Atmospheric 
Distributed  Data  System  (ADDS),  being  currently  tested  on  Air  Force  air¬ 
craft,  is  necessary.  With  improvements  in  communications,  navigation, 
data  density,  and  data  accuracy,  the  analysis  and  predictive  models 
surely  will  improve. 

The  results  of  this  study  do  not  satisfy  totally  the  original  ob¬ 
jectives  as  stated  when  the  research  was  started  (see  Chapter  II);  how¬ 
ever,  using  a  more  reliable  data  base  as  made  available  by  NOAA  research 
aircraft  or  the  future  ADDS  program  may  show  significant  'mprovement. 

Of  particular  interest  will  be  the  improvement  of  Vg^  computations  with 
greater  navigational  precision  and  data  density.  Even  with  the  likeli¬ 
hood  of  continued  wind  measurement  errors  within  the  tropical  cyclone 
environment,  significant  gains  should  be  possible  in  the  near  future  if 
forward  steps  are  taken.  This  program  and  concept  of  analysis  is  offer¬ 
ed  for  future  study  using  improved  data  bases  soon  to  be  available. 
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APPENDIX 


This  Appendix  discusses  the  data  handling  and  processing  procedures 
used  in  determining  the  intensity,  strength,  and  size  of  the  tropical 
storms  evaluated.  Information  on  the  source  of  the  original  data  is 
presented.  The  format  in  which  the  data  are  ordered  for  computation  is 
described.  The  equations  used  in  calculating  are  derived  and  explained, 
including  descriptions  of  how  many  of  the  quantities  are  determined. 
Finally,  the  detailed  computer  program,  including  annotations,  is  dis¬ 
played. 

The  data  comprising  the  original  aircraft  reconnaissance  observa¬ 
tions  were  recorded  on  standard  forms  and  archived  in  the  National 
Climatic  Data  Center,  Asheville,  North  Carolina.  The  forms  for  1980- 
1981  were  obtained  from  the  U.S.  Air  Force  Environmental  Technical 
Applications  Center.  The  data  from  the  forms  were  put  in  a  coded  format 
which  could  be  used  in  the  computer  process.  Because  of  the  time  in¬ 
volved  in  this  phase,  only  one  year  of  data  (1980)  and  the  data  for 
Super  Typhoon  Tip  (1979)  were  used  for  this  study. 

The  data  were  coded  in  the  following  format: 

aa  bbccc  ddd  eeee  ffggg  hhhh  iijjj  kkkk. 

The  significance  of  each  group  is: 

aa  -  code  of  observation  level  and  type 

bb  -  consecutive  day  after  first  mission  flown  on  storm 

ccc  -  Greenwich  Mean  Time  (GMT)  of  observation 

ddd  -  latitude  of  observation  in  degrees  and  tenths 

eeee  -  longitude  of  observation  in  degrees  and  tenths 

ff  -  flight  level  wind  direction  n  tens  of  degrees 

ggg  -  flight  level  wind  speed  in  knots 

hhhh  -  height  of  standard  surface  in  meters 


ii  -  surface  wind  direction  in  tens  of  degrees 
jjj  -  surface  wind  speed  in  knots 
kkkk  -  sea  level  pressure  in  millibars 

All  data  were  ordered  in  integer  format  for  simplicity  and  effi¬ 
ciency  of  computation.  Further  explanation  of  the  coded  data  is  inclu¬ 
ded  in  the  commentary  portion  of  the  Appendix. 

Initially,  data  were  taken  from  the  input  data  base  which  was  stor¬ 
ed  on  disk  in  the  format  described  in  Chapter  IV  and  in  the  Appendix. 
Three  separate  output  data  sets  were  generated  and  labeled  according  to 
storm  name,  mission  number,  and  number  of  observations  for  the  mission. 
Observations  for  an  entire  mission  were  placed  into  memory  at  one  time. 
The  times  of  these  observations  were  changed  from  the  GMT  format  to  a 
continuous  minute  format  starting  at  0000  GMT  of  the  first  day  that  a 
mission  was  flown  on  a  particular  storm.  This  helped  reduce  the  comput¬ 
ing  time  required  for  determinations  of  system  velocity  and  center 
position. 

As  the  time  was  converted,  sea  level  pressures  were  extrapolated 
for  observations  taken  at  the  70  kPa  level.  This  allowed  both  a  sea 
level  pressure  gradient  and  a  70  kPa  height  gradient  to  be  computed. 
Jordan's  (1958)  mean  summer  tropical  atmosphere  was  adjusted  for  temper¬ 
ature  and  relative  humidity  differences  for  the  innermost  220  km  of 
western  Pacific  typhoons  (Table  A-1)  after  Gray  ^  al_.  (1975).  This 
modified  atmosphere  (Fig.  A-1)  gives  a  standard  100-70  kPa  thickness  of 
3063  m  for  the  storm  environment.  Using  this  value  plus  the  given  70 
kPa  height  (H^q)  for  each  observation  outside  of  the  center,  a  100  kPa 
height  value  (Hiqq)  was  computed. 
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Table  A-1 .  Mean  summer  tropical  atmosphere  of  Jordan  (1958) 
adjusted  for  the  western  Pacific  typhoon  environment  (from  Gray  et  al . , 
1975). 


Level  (kPa)  T  (°C)  RH  {%)  (°C)  w  (g/kg)  (°C) 


^Tv  of  adjusted  mean 
tropical  atmosphere 
aT  of^  the  adjusted  mean 
^tropical  atmosphere 

4  T  of  the  standard 
a  atmosphere _ 


r  o: 

'0> 


10  20  30  A0°C 

TEMPERATURE 

A  tephigram  which  compares  ambient  and  virtual  temperature 
profiles  of  tropical  and  standard  atmospheres.  These  profiles 
include  the  virtual  temperature  (•  )  and  ambient  temperature 
(^)  profiles  of  the  adjusted  mean  tropical  atmosphere  of  the 
western  Pacific  as  well  as  the  ambient  temperature  profile  of 
the  standard  atmosphere  C^)  (after  Gray  et  al . ,  1975). 


The  hypsometric  equation  was  used  in  the  form 


where  C  is  the  hypsometric  constant  and  p^  is  the  sea  level  pressure  in 
kPa.  The  hypsometric  constant  is 

C  =  g/Rjj  ,  (A-3) 

where  g  is  the  gravitational  constant,  is  the  gas  constant  for  dry 
air,  and  is  the  assumed  mean  virtual  temperature  between  the  100  kPa 
level  and  the  surface.  The  modified  atmosphere  of  Table  A-1  indicates 
that  T^  is  29.1°C  at  100  kPa  and  this  is  a  realistic  approximation  of  T^ 
for  most  cases.  Using  this  to  establish  C  and  rearranging  givrs 

Pq  =  100  /  [exp 

as. the  sea  level  pressure  in  the  tropical  cyclone  environment  outside  of 
the  center. 

When  a  fix  was  made  of  the  70  kPa  level  circulation  center  without 
the  benefit  of  a  minimum  sea  level  pressure  observation  at  the  surface 
center  (p^.).  minimum  sea  level  pressure  was  extrapolated  from  70  kPa 
using  the  70  kPa  center  minimum  height  and  the  Jordan  (1957)  formula. 

Pc  =  .115(H7q)  +  645  .  (A-4) 

This  relationship  works  well  even  when  surface  and  70  kPa  centers  are 
separated  horizontally,  even  a  significant  distance.  It  is  preferred  to 
the  former  method  of  computation  because  of  the  variability  of 
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temperature  and  100-70  kPa  thicknesses  in  the  center  environment,  where¬ 
as  these  elements  are  fairly  constant  outside  of  the  center  of  most 
storms . 

Since  the  program  may  be  used  in  near  real  time  situations,  a  means 
of  determining  system  velocity  and  position  needed  to  be  developed  using 
only  reconnaissance  data,  without  relying  on  the  best  track  or  forecast 
values.  The  goal  was  to  make  the  program  as  self-sufficient  as  pos¬ 
sible.  In  the  second  part  of  the  program,  the  data  were  searched  for 
flight  level  and  surface  center  observations.  The  times  and  positions 
of  these  observations  were  then  used  to  compute  the  movement  of  the 
system.  The  velocity  of  the  system  was  computed  by  determining  the 
changes  in  latitude  and  longitude  between  fixes.  With  the  fix  posi¬ 
tions,  fix  times,  and  system  velocities  in  computer  memory,  a  center 
position  could  be  computed  for  any  time. 

In  part  3  of  the  program,  the  following  computations  were  neces¬ 
sary.  First,  a  determination  was  made  of  the  difference  between  the 
current  observation  time  and  the  next  future  fix  time  (or  the  last  fix 
time  of  the  storm).  Once  this  was  determined,  the  system  fix  positions 
and  velocities  that  were  stored  in  computer  memory  were  used  to  compute 
the  current  center  positions  (both  flight  level  and  surface).  Each  in¬ 
dividual  observation  was  then  plotted  relative  to  the  current  center 
position. 

In  order  to  compute  wind  and  pressure  gradient  relationships,  ob¬ 
servations  must  be  normalized  to  the  pressure  field.  If  assumption  is 
made  that  (for  the  tropical  cyclone  environment)  the  pressure  field  is 
symmetrical  about  the  center,  then  the  normalization  can  be  accomplished 


by  transposing  observations  to  a  common  radial  line  from  the  center. 

The  assumption  of  pressure  symmetry  is  realistic  for  the  tropical 
cyclone  mission,  considering  the  small  distances  between  most  observa¬ 
tions. 

.  In  Fig.  A-2,  a  situation  where  two  consecutive  observations  undergo 
transposition  is  seen.  Observation  1  is  positioned  at  angle  and 
radial  distance  r-j  at  time  t-j .  Observation  2  is  positioned  at  angle  02 
and  radial  distance  r2  at  time  t2.  It  can  be  seen  that  these  two  obser¬ 
vations  do  not  initially  fall  on  a  common  radial  line  from  the  center 
unless  we  transpose  observation  1  to  a  position  at  distance  r^  along  the 
radial  line  through  observation  2  and  the  center  of  the  cyclone.  By 
rotating  observation  1  through  an  angle,  p  =  02  -  0i ,  it  arrives  along 
the  radial  line  through  observation  2  at  radial  distance  r.|  and  time  t2. 
During  this  process,  the  pressure  of  observation  1  (p^)  remains  con¬ 
stant  since  we  assumed  symmetry.  However,  the  wind  direction  of  the 
transposed  observation  must  be  adjusted  for  the  angle  of  transposition 
(e).  Now  a  realistic  pressure  gradient  can  be  attained  since  ap  is 
normal  to  the  symmetric  pressure  field. 

Observed  winds  must  also  be  corrected  for  system  motion  before  be¬ 
ing  used  in  the  computations.  The  absolute  velocity  of  the  wind  con¬ 
sists  of  that  portion  due  to  the  pressure  gradient  and  the  component  due 
to  the  system's  motion.  To  adjust  the  winds  to  a  relative  velocity,  a 
sub-program  was  written  which  split  each  observed  wind  velocity  into  u 
and  V  components  and  subtracted  the  u  and  v  component  of  the  system 
motion.  The  transposition  of  observation  1  was  then  made  through  angle 
6  with  a  corresponding  correction  to  the  wind  direction. 


Fig.  A-2.  The  transposition  of  observations  to  a  common  radial  line 
normal  to  a  symmetrical  pressure  field.  Observation  1  at 
angle  cfi,  radial  distance  ri,  and  pressure  Pi,  is  transposed 
to  the  same  radial  line  as  observation  2  at  angle  ({)2,  radial 
distance  r2.  and  pressure  Pj.  The  angle  S  is  the  angle  of 
transposition  of  observation  1  to  the  same  radial  line  as 
observation  2. 


65 


When  consecutive  observations  contained  pressure/height  data,  the 
pressure/height  gradient  was  compared  to  the  observed  winds.  Now  that 
the  observations  and  winds  are  transposed  to  a  common  radial,  the  gradi¬ 
ent  and  winds  may  be  correlated.  To  do  this,  the  maximum  wind  between 
the  two  observations  as  well  as  the  average  wind  were  computed.  These 
data  were  then  interpolated  to  a  midpoint  position  between  the  two  ob¬ 
servations;  they  represent  the  average  and  maximum  winds  which  are 
associated  with  the  pressure/height  gradient  at  radial  distance  r^ 

(Fig.  A-2). 

On  occasion,  consecutive  observations  do  not  contain  pressure/ 
height  data.  This  occurs  most  frequently  on  low  level  missions  where 
midpoint  winds  are  taken  between  pressure  observations.  When  this 
occurred,  the  program  transposed  the  previous  observations  with 
pressure/height  data  and  one  intermediate  wind  observation  (Fig.  A-3). 
All  three  transposed  winds  were  evaluated  to  determine  average  and  maxi¬ 
mum  winds,  which  in  turn  were  correlated  to  the  gradient. 

It  must  be  remembered  that  simultaneous  correlations  of  flight 
level  and  surface  level  winds  and  gradients  were  computed.  This  was 
necessary  because  of  the  frequent  vertical  tilt  observed  in  tropical 
cyclones  (Fig.  A-4).  Problems  arise  when  we  extrapolate  sea  level 
pressures  from  the  70  kPa  level  and  try  to  infer  gradient  relationships 
under  these  circumstances.  It  is  safe  to  assume  that  under  most  circum¬ 
stances  when  vertical  tilt  is  not  extreme,  the  sea  level  pressure  extra¬ 
polation  is  fairly  accurate  and  the  sea  level  pressure  gradient  (SLPG) 
can  be  computed  with  reasonable  accuracy.  However,  for  locations  within 
a  one  degree  of  latitude  radius  of  the  vortex  center,  the  tilt  of  the 
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n  Fig.  A-2,  except  with  an  intermediate  observation 
ervation  2)  without  pressure  data.  The  angle  6'  is 
angle  of  transposition  of  the  intermediate  observation 
he  same  radial  line  as  observation  3. 


UPPER  VORTEX 


Fig.  A-4.  Vertical  tilt  of  a  tropical  cyclone  vortex  between  two 
pressure  levels.  The  observation  point  (X)  is  a  certain 
distance  (r^J)  from  the  vortex  center  at  the  flight  pressure 
level  (p^j)  and  a  greater  distance  (r^)  from  the  surface 
vortex  center.  The  surface  center's  pressure  level  (p^) 
is  shown  as  it  slopes  upward  with  distance  from  the  center. 


vortex  may  significantly  influence  the  actual  SLPG.  When  the  tilt  is 
extreme,  the  accurate  determination  of  SLPG  is  further  complicated,  yet 
under  these  circumstances,  the  tropical  cyclone  is  usually  poorly 
organized,  and  the  winds  are  generally  light. 

With  maximum  and  average  flight  level  winds  available  for  nearly 
every  set  of  observations,  correlations  between  these  winds  and  the 
flight  level  height  gradient  (FLHG)  were  determined.  If  surface  winds 
were  obsc  vable,  correlations  were  drawn  between  the  FLHG  and  the  maxi¬ 
mum  and  average  surface  wind.  In  addition,  with  sea  level  pressure 
extrapolated  from  70  kPa  heights,  surface  wind  correlations  with  the 
SLPG  were  possible. 

After  these  correlations  were  determined,  data  were  distributed  to 
two  of  the  three  output  files.  Output  file  number  one  contained  the 
input  data  with;  1)  the  time  converted  to  continuous  minute  format, 

2)  the  latitude  and  longitude,  3)  the  flight  level  and  surface  winds 
corrected  for  system  motion,  and  4)  the  sea  level  pressure  (extrapolated 
from  70  kPa  heights  if  necessary).  The  second  output  file  contained 
distances  from  the  flight  level  and  surface  centers  to  the  midpoint  be¬ 
tween  observations  with  pressure/height  data.  It  also  listed  the  total 
correlations,  which  consisted  of:  1)  average  flight  level  wind  with  the 
FLHG,  2)  average  flight  level  wind  with  the  SLPG,  3)  maximum  observed 
flight  level  wind  with  the  FLHG,  4)  maximum  observed  flight  level  wind 
with  the  SLPG,  5)  average  surface  wind  with  the  FLHG,  6)  maximum  observ¬ 
ed  surface  wind  with  the  FLHG,  7)  average  surface  wind  with  the  SLPG, 
and/or  8)  maximum  observed  surface  wind  with  the  SLPG. 

If  the  results  prove  meaningful,  the  analysis  of  tropical  cyclone 


intensity  and  strength  characteristics  may  be  done  using  the  wind  and 
gradient  correlations  described  above.  Intensity  may  be  approximated  by 
determining  the  maximum  gradient  and  relating  it  to  a  maximum  wind 
value.  Strength  may  be  determined  by  averaging  the  gradient  at  speci¬ 
fied  radial  distances  from  the  cyclone's  center.  The  averaging  would  go 
out  to  the  limit  of  the  cyclonic  circulation,  and  an  average  strength 
could  be  computed.  However,  in  order  to  analyze  tropical  cyclone  size, 
the  program  must  go  further. 

First,  a  definition  of  size  must  be  determined.  Ideally,  the  size 
is  determined  by  the  average  radius  outside  which  any  tropical  cyclone 
effect  is  negligible.  Practically,  this  has  been  done  in  the  past  by 
evaluating  the  average  radius  of  the  outermost  closed  isobar  (Merrill, 
1982).  An  analysis  of  this  type  may  be  done  provided  that  the  synoptic 
and  ship  reports  will  allow  an  accurate  pressure  analysis.  However,  in 
some  parts  of  the  ocean,  this  is  hard  to  accomplish. 

One  of  the  objectives  is  to  determine  tropical  cyclone  size  in  near 
real  time.  The  best  approach  incorporates  the  Holland  (1980)  hurricane/ 
typhoon  wind  and  pressure  profile  equation. 


•  Pc^  /  (P  ■  ^  ’ 


(A-5) 


where  r  is  the  radial  distance  from  the  center,  p^  is  the  environmental 
pressure  (theoretically  at  infinite  radius),  p  is  the  central  pressure, 
p  is  the  pressure  at  radius  r,  and  A  and  B  are  scaling  parameters.  Re¬ 
arranging  this  equation  yields 


r  =  {A  /  In  [(p„  -  p^)  /  (p  -  p^)]}' 


(A-6) 


Fig,  A-5  approximates  the  family  of  rectangular  hyperbolas  which  repre¬ 
sent  the  pressure  profiles  of  various  typhoons.  If  p'  is  substituted 
for  p  in  (A-6)  and  set  at  the  tropical  cyclone's  assumed  outermost 
closed  isobar  value,  then  we  can  substitute  r'  for  r  as  the  storm's 

radius  and  solve  for  r'  once  A  and  B  are  established.  The  equation  then 

becomes 

r'  =  {A  /  In  [(p^  -  p^)  /  (p-  -  .  (A-7) 

The  scaling  parameter  B  determines  the  shape  of  the  profile  and 
varies  between  about  1.0  and  2.5  for  most  typhoons.  The  stronger  the 

pressure  gradient  at  the  radius  of  maximum  wind,  the  higher  the  value  of 

B.  For  this  study,  Holland's  minimum  limit  of  B  has  been  lowered  from 
1.0  to  0.5  in  order  to  consider  the  weaker  pressure  gradients  of  sub¬ 
typhoon  systems.  In  the  program,  B  can  best  be  approximated  by  measur¬ 
ing  the  innermost  (0-55  km)  pressure  gradient  (IPG)  and  translating  it 
into  a  corresponding  value  of  B.  Using  the  observed  extremes  (weak  and 
strong)  of  IPG  for  the  1980  data,  a  sub-program  was  developed  which 
could  interpolate  a  value  of  B  from  between  0.5  and  2.5  for  every  fix 
(measured  IPG). 

Once  B  is  determined,  a  value  of  A  may  be  determined  by  solution  of 
(A-5)  with  values  of  r  and  p  known.  Solving  for  r'  in  (A-7)  then  gives 
the  storm's  radial  size. 

Representative  values  of  p^  and  p'  are  critical  to  the  accurate 
computation  of  size  (r').  A  method  to  determine  objectively  p^  and  p' 
was  included  in  the  program.  However,  as  discussed  later,  this  may  be 
done  best,  subjectively.  The  objective  method  makes  use  of  the  Atlas  of 


RADIAL  DISTANCE  FROM  THE  CENTER  (km) 


Pressure  profiles  of  randomly  selected  typhoons  used  in 
this  study.  The  parameter  (p  -  Pc)/(Pn  -  Pc)  is  used 
to  normalize  the  profiles  for  variations  due  to  different 
central  and  environmental  pressures. 
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Mean  Sea-Level  Pressure  (Chary,  1982)  for  the  determination  of  mean 
environmental  pressures  at  latitudes  equatorward  of  30°,  Environmental 
pressure,  as  used  here,  is  defined  as  the  pressure  at  infinite  radius 
from  the  storm.  Fig.  A-6  shows  a  typical  SLP  pattern  for  the  western 
North  Pacific  in  October  (peak  of  the  tropical  cyclone  season).  As 
would  be  expected  from  the  theory  of  the  general  circulation  of  the 
atmosphere,  low  pressure  near  the  equator  increases  to  high  pressure 
near  30  degrees  latitude  (subtropical  ridge)  and  then  decreases  toward 
60  degrees  latitude.  A  value  of  101  kPa  is  a  reasonable  year  round 
estimate  of  environmental  pressure  in  the  equatorial  region.  However, 
at  30°N  the  pressure  varies  seasonally  as  well  as  longitudinally.  We 
can  suppose  that  the  environmental  pressure  at  30°N  is  at  least  as  high 
as  the  central  pressure  of  the  subtropical  ridge.  With  a  tropical 
cyclone  breaking  through  the  ridge  axis,  the  east-west  pressure  profile 
will  be  largely  determined  by  the  central  strength  of  the  subtropical 
ridge.  Additionally,  in  the  mid-winter  months,  the  strenghtening  of  the 
continental  high  pressure  of  Asia  may  increase  environmental  pressures 
at  30°N  to  over  102.7  kPa.  Because  of  this  and  other  factors  (to  be 
described  later),  a  value  of  102.7  kPa  was  chosen  as  the  environmental 
pressure  at  30°N.  The  value  of  p^  can  then  be  expressed  as  a  sinusoidal 
function  of  latitude  with  an  amplitude  of  1.7  kPa  between  the  equator 
and  30°N. 

To  demonstrate  a  point,  two  hypothetical  cyclones  are  positioned  as 
shown  in  Fig.  A-6.  One  cyclone  (A)  is  breaking  through  the  subtropical 
ridge  on  a  northerly  track,  and  the  other  cyclone  (B)  remains  embedded 
in  the  easterlies  south  of  the  subtropical  ridge.  Evaluating  an  east- 


Fig.  A-6.  Two  tropical  cyclones  positioned  in  a  typical  western  North 
Pacific  sea  level  pressure  field  for  October  (after  Chary, 
1982).  Isobars  are  labeled  in  kPa. 


west  cross  section  of  pressure  for  cyclone  A  shows  that  increases  to 
over  102.4  kPa  in  China  and  to  over  102.2  kPa  in  the  central  Pacific.  A 
north-south  cross  section  of  the  same  cyclone,  on  the  other  hand,  shows 
a  significantly  smaller  increase  (to  about  101.2  kPa  in  the  north  and  to 
101.1  kPa  in  the  south).  Using  the  sinusoidal  relationship  discussed 
previously,  the  computer  program  determines  p^  to  be  102.2  kPa  (based  on 
the  latitude  of  the  cyclone),  which  would  compare  well  with  the  east- 
west  cross  section.  Evaluating  cyclone  B  in  the  east-west  cross  section 
indicates  that  p^  increases  to  over  101.9  kPa  in  the  central  Pacific  and 
to  over  101.6  kPa  to  the  west.  Evaluating  the  north-south  cross  section 
shows  a  p^  of  nearly  102.1  kPa  to  the  north  and  only  101.0  kPa  to  the 
south.  The  computer  program  determines  that,  for  the  latitude  of 
cyclone  B,  p^  is  102.1  kPa,  which  may  be  too  high  for  all  but  the  north¬ 
ern  quadrant  in  this  case,  but  ft  is  climatologically  reasonable. 

Although  subjective  determination  of  p^  may  be  the  best  procedure 
in  some  cases,  an  adequate  subjective  analysis  of  the  environmental, 
pressure  is  not  always  available  in  real  time.  Added  to  this,  as  will 
be  discussed  later,  is  the  fact  that  maintaining  a  steep  rate  of  in¬ 
crease  between  p^  and  p^  is  advantageous  in  size  computations,  and  the 
program  determinations  of  p^  tend  to  favor  these  steep  increases. 

As  program  execution  was  running,  no  size  computations  were  per¬ 
formed  until  the  initial  center  fix  for  any  particular  mission  was  made. 
At  that  time,  a  measurement  of  the  IPG  was  available  for  the  approxima¬ 
tion  of  B.  During  every  fix  thereafter,  a  new  measurement  of  IPG  and 
approximation  of  B  was  made  until  the  end  of  the  mission.  Once  the 
initial  approximation  of  B  was  available,  computations  of  r'  were  made 


after  every  observation  where  pressure  gradient  was  obtained  as  long  as 
the  radial  distance  from  the  storm's  center  to  the  observation  was 
greater  than  that  of  a  previous  computation.  To  do  this,  the  scaling 
parameter  A  was  solved  in  (A-5)  using  the  observation's  radial  distance 
as  r,  the  observation's  sea  level  pressure  as  p,  and  a  time  interpolated 
or  extrapolated  value  of  central  sea  level  pressure  based  on  the  known 
fix  values  for  p^.  Equation  (A-7)  was  then  solved  for  r'  with  p'  set  at 
an  assumed  value. 

The  ideal  method  of  determining  tropical  cyclone  size  is  to  measure 
it  directly.  As  this  cannot  be  done  in  real  time  for  most  storms,  the 
use  of  (A-7)  was  the  best  method  of  approximation.  Fig.  A-7  shows  that 
if  the  observation  point  is  outside  the  radial  limit  (or  outermost  clos¬ 
ed  isobar)  of  the  storm  (i.e.,  r  >  r')  and  if  all  data  indicate  that  p' 
is  set  accurately,  then  interpolation  can  be  used  to  determine  accurate¬ 
ly  where  p'  is  (in  terms  of  r')  along  a  radial  from  the  center.  How¬ 
ever,  if  p'  is  not  accurate,  r'  will  be  inaccurate  as  well  due  to  the 

shallow  slope  of  the  curve.  If  r  <  r',  then  extrapolation  must  be  used 

to  estimate  the  size.  This  can  be  very  difficult,  because  very  small 

errors  in  p^  and/or  p'  can  cause  very  large  errors  when  r'  is  extra¬ 
polated,  unless  r  is  close  enough  to  r'  so  that  the  extrapolation  is 
over  only  a  very  short  distance.  For  these  reasons,  only  the  computa¬ 
tion  made  at  the  greatest  distance  from  the  storm's  center  was  used  in 
size  computations. 

As  can  be  seen  in  Fig.  A-7,  if  p'  was  incorrectly  computed  too 
large  and  the  outermost  observed  pressure  (p)  was  within  p'  (r  <  r'), 
then  r'  may  be  grossly  inflated.  Effort  has  been  made  to  avoid  this 


observations  used  to  define  the  outer  limit  pressure  (p')  and  size 

(r')  of  a  tropical  cyclone.  In  (a),  the  radial  distance  (r')  of  tlie 

assumed  outeriiicst  closed  isobar  value  (p')  is  interpolated  along  the 
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circumstance  by  carefully  adjusting  p'  in  the  program.  The  program 
started  by  assuming  p'  was  100.8  kPa.  This  is  a  reasonable  value  to  ex¬ 
pect  for  the  outermost  closed  isobar  on  synoptic  analyses  of  western 
North  Pacific  tropical  cyclones.  Using  the  combined  1980-1981  data  from 
non-storm  related  observations  of  flight  level  and  surface  wind,  mean 
characteristic  values  for  the  western  North  Pacific  are  13  and  11  m  s"^ 
respectively.  If  both  flight  level  and  surface  winds  dropped  below 
these  threshold  values  at  r  (with  pressure  p)  when  p  <  100.8  kPa,  it  was 
assumed  that  p'  =  p.  Conversely,  if  p  >  100.8  kPa  and  if  either  flight 
level  or  surface  winds  were  above  18  or  13  m  s~^  respectively,  then  the 
winds  were  considered  to  be  influenced  by  the  cyclone,  p'  =  p,  and  the 
process  continued  for  larger  values  of  p  and  r  until  flight  level  and 
surface  winds  dropped  below  the  18  and  13  m  s~^  thresholds.  Since  the 
computed  wind  velocities  in  the  program  had  system  velocity  subtracted, 
velocities  above  or  below  the  thresholds  were  used  to  indicate  the 
presence  or  absence  of  tropical  cyclone  effect  and  to  adjust  p'. 

In  addition  to  the  above  procedure,  p^  was  never  allowed  to  be 
closer  than  0.3  kPa  from  p'.  This  helped  avoid  extremely  shallow  slopes 
between  p'  and  p^  at  low  latitudes  and  helped  prevent  inflation  problems 
when  extrapolating  r'. 

Even  when  precautions  were  taken  to  avoid  inflation  of  r'  when 
extrapolating,  the  majority  of  all  tropical  cyclone  missions  were  flown 
totally  within  r',  and  therefore,  extrapolation,  and  consequently  infla¬ 
tion,  does  occur.  The  effects  of  this  inflation  are  most  noticeable 
when  successive  size  computations  are  observed  to  exhibit  a  large 
variance.  In  order  to  avoid  this  variation,  the  reconnaissance  data  may 


be  supplemented  by  a  ship  report  (if  available)  at  some  distance  (rg) 
greater  than  the  outermost  reconnaissance  observation  and  preferably 
greater  than  r'.  Provided  that  the  observation  time  is  reasonably 
coincident  with  the  reconnaissance  observations  in  the  same  quadrant  of 
the  storm,  the  observation  may  be  added  in  with  the  data  and  will  allow 
r'  to  be  interpolated  instead  of  extrapolated. 

Subjective  substitution  of  p'  and  p^  will  help  smooth  the  ir¬ 
regularities  of  this  rather  mechanical  program  size  computation.  The 
program  has  a  tendency  to  exaggerate  the  impact  of  the  wind  thresholds 
described  earlier  in  determining  p'.  For  example,  considering  tropical 
cyclone  B  in  Fig.  A-6,  the  resulting  pressure  field  is  highly  asymmetric 
with  stronger  gradients  and  winds  in  the  northern  quadrant.  If  the  only 
track  into  and  out  of  the  storm  is  from  the  south,  even  a  storm  with 
20  m  s'^  winds  to  the  north  may  possess  less  than  10  m  s”^  winds  well  in 
toward  the  center  on  the  south  side,  even  though  the  actual,  average 
storm  radius  extends  outward  a  considerable  distance.  Under  these 
circumstances  the  program  tended  to  decrease  p'  toward  p^  until  11  m  s”^ 
surface  winds  or  13  m  s”^  flight  level  winds  were  encountered.  Without 
the  counter  effect  of  a  track  to  the  north  of  the  storm,  the  computed 
radius  of  the  storm  would  be  an  underestimate. 

The  calculations  and  procedures  described  in  the  preceding  pages 
were  combined  into  a  computer  program  developed  and  used  to  determine 
tropical  cyclone  vortex  motion  and  vertical  tilt  as  well  as  cyclone 
strength,  intensity  and  size.  The  program  also  allowed  computation  of 
desired  pressure/wind  relationships.  The  program  was  written  in  basic 
Fortran  for  the  Harris  100  computer  at  the  Department  of  Meteorology, 
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AUTHOR:  CHARLES  B.  STANFIELD 

PROGRAM;  TROPICAL  CYCLONE  WIND/PRESSURE  STUDY 
PURPOSE:  THIS  PROGRAM  PLOTS  WINDS  AND  PRESSURES  VERSUS  RADIAL 
DISTANCE  FROM  THE  CENTER  OF  A  TROPICAL  CYCLONE  VORTEX.  IT  WILL 
DETERMINE  CHARACTERISTIC  PRESSURE/HEIGHT  GRADIENTS  AND  RELATE 
THEM  TO  THE  OBSERVED  WIND  FIELD.  THESE  RELATIONSHIPS  CAN  THEN  BE 
USED  TO  QUANTIFY  TROPICAL  CYCLONE  INTENSITY,  STRENGTH.  AND  SIZE. 

A  BRIEF  DESCRIPTION  WILL  PRECEDE  EACH  PART  OF  THIS  PROGRAM.  AND  A 
VARIABLE  LIST  WILL  ALSO  BE  GIVEN. 


INTEGER/REAL/CHARACTER/DATA  STATEMENTS* 


INTEGER  CQDEOO.SO)  .TIME (30,50)  ,FLW( 30.50 )  .HSSOO,  50)  .SFW(  30. 50)  , 
.  SLP(30.50),SFLG.SM1N1.LCMIN(60).SCMIN(60),DIFTIM,N(30).NUM(30), 

.  F L ANG 1 . F L ANG2 . F LANG3 , SLANG1 , SLANG2 , SL ANG3 , SCSLP ( 60 ) 

REAL  LAT(30.50) .L0N(30.50).LLAT1 .LL0N1 .LVELAT,LVEL0N,LCLAT(60) , 

.  LCL0N(60) , LCVLAT(60) .LCVL0N(60) .SCLAT(60) .SCL0N(6O) .SCVLAT(60) . 

.  SCVL0N(60) ,DISLLA(50) .DISLL0(5O).DISSLA(5O).DISSL0(5O).LCKLAT. 

.  LCKLON 

character*?  NAME 

DATA  K1 ,K2.K3.C1 , C2 . C3 . C4 . C5/ 1000 . 10.90000,57.296.3598. 14. 
.59.969.3063. , .000113/ 


PART  1 

THE  DATA  ARE  FILED  UNDER  INPUT.  THE  DATA  BEGINS  WITH  A  CODE  TO 
DESIGNATE  THE  STORM  (I.E.  B0000010NE  6).  THIS  EXAMPLE  DESIG¬ 

NATES  THE  STORM  TO  BE  FROM  1980  (80).  WITH  NUMERICAL  DESIGNATION 
ONE  (00001)  AND  NAME  "ONE".  SIX  SEPARATE  RECONNAISSANCE  DATA  SETS 
FOLLOW  (6).  EACH  SEPARATE  DATA  SET  BEGINS  WITH  A  CODE  (I.E. 
8002001  20),  WHERE  THE  "02"  REPRESENTS  THE  MISSION  NUMBER. 

THE  "20"  REPRESENTS  THE  NUMBER  OF  OBSERVATIONS  ON  THAT  MISSION. 
THE  FIRST  GROUP  OF  NUMBERS  ON  THE  DATA  LINE  FOR  THE  INDIVIDUAL 
OBSERVATION  IS  A  CODE  FOR  THE  TYPE  OF  OBSERVATION  (I.E.  10*700MB 
OBSERVATION  TAKEN  OUTSIDE  THE  CENTER,  11=L0W  LEVEL  (1500FT)  OB¬ 
SERVATION  TAKEN  OUTSIDE  THE  CENTER,  20'SURFACE  CENTER  OBSERVATION 
TAKEN  FROM  700MB.  21*SURFACE  CENTER  OBSERVATION  TAKEN  FROM 
1500FT.  30«700MB  CENTER  OBSERVATION  TAKEN  AT  700MB,  AND  41»L0W 
LEVEL  (1500FT)  CENTER  POSITION  TAKEN  AT  1500FT).  THE  SECOND  GROUP 
OF  NUMBERS  INDICATES  THE  TIME  (I.E,  010035:  DAY  ONE  (01)  OF  DATA 
gathered  on  the  storm  at  time  0035  GMT).  THE  THIRD  AND  FOURTH 
GROUPS  REPRESENT  THE  LATITUDE  AND  LONGITUDE  IN  TENTHS  OF  DEGREES 
RESPECTIVELY,  THE  FIFTH  GROUP  IS  WIND  DIRECTION  AND  SPEED  OF  THE 
flight  level  wind  (I.E  17017;  170  DEGREES  AT  17  KNOTS).  THE 
SIXTH  GROUP  IS  The  height  of  the  700MB  SURFACE  (WHEN  GIVEN).  THE 
seventh  group  is  the  ESTIMATED  SURFACE  WIND  AND  THE  FINAL  GROUP 
IS  SEA  level  pressure  (SLP)  in  millibars 


THE  FIRST  PART  OF  THE  PROGRAM  WILL  READ  THE  INDIVIDUAL  STORM 
number,  name,  and  the  number  of  reconnaissance  MISSIONS  FLOWN,  IT 
WI>.  READ  AND  P..ACE  INTO  MEMORY  THE  INDIVIDUAL  MISSION  NUMBER,  THE 
NUMBER  Of  OBSERVATIONS  FOR  EACH  MISSION.  AND  THE  INDIVIDUAL  OBSER- 
VA’-;orYf  fQR  the  entire  STORM  THE  CODED  TIME  WILL  BE  CONVERTED 
INTC  ..  CONSECUTIVE  MINUTE  STARTING  AT  OOOOZ  ON  THE  FIRST  DAY  THAT 
data  was  gathered  ON  the  storm,  this  time  can  THEN  BE  USED  IN 
SYSTEM  velocity  computations,  the  data  is  also  manipulated  to 
extrapolate  SLP  FROM  700MB  DATA.  IF  NECESSARY. 


VARIABLE  list- 


KLL  -  flag  for  1500FT  FLIGHT  LEVEL 
KHH  -  FLAG  FOR  700MB  FLIGHT  LEVEL 
NU  -  STORM  NUMBER  ID 
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NUM(K)  -  MISSION  NUMBER  ID 

NAME  -  STORM  NAME 

NN  -  NUMBER  OF  MISSIONS  ON  STORM 

N(K)  -  NUMBER  OF  OBS  IN  MISSION 

K  -  DO  LOOP  COUNTER  FOR  MISSION  NUMBER 

I  -  DO  LOOP  COUNTER  FOR  OBS  NUMBER 

C0DE(K,I)  -  OBS  LEVEL  CODE  FOR  OBS  I  OF  MISSION  K 

TIME(K.I)  -  TIME  OF  OBS  I 

LAT(K,I)  -  LATITUDE  OF  OBS  1 

LON(K.I)  -  LONGITUDE  OF  OBS  I 

FLW(K,I)  -  FLIGHT  LEVEL  WIND  OF  OBS  I 

HSS(t<,I)  -  HEIGHT  OF  STANDARD  SURFACE  FOR  OBS  I 

SFW(K,I)  -  SURFACE  WIND  FOR  OBS  I 

SLP(K,I)  -  SEA  LEVEL  PRESSURE  FOR  DBS  I 

ITIM  -  TEMPORARY  TIME 

MIN  -  TEMPORARY  MINUTE  OF  OBS  FROM  DAY  ONE  OF  THE  STORM 
J  -  COUNTER  OF  THE  FIXES 
DZ  -  HEIGHT  OF  1000MB  SURFACE 

C4  -  STANDARD  1000/700MB  THICKNESS  FDR  STORM  ENVIRONMENT 

25  -  hypsometric  CONSTANT  USING  VIRTUAL  TEMPERATURT  CONSTANT  29. 1C 


98  READ( 16, 1 , END=99)NU.NAME.NN 

1  F0RMAT(I7.A7, 12) 

WRITE( 17,2)NU.NAME,NN 
WRITE) 18.2)NU.NAME.NN 
WRITE) 19,2)NU,NAME.NN 

2  FORMAT)  IX , 17 , A7 , 12  ) 

KLL  =  0 

KHH  =  0 
J  =  0 

D07  K»1 ,NN 

READ) ie.9)NUM)K),N)K) 

9  FORMAT)I7,7X,I2) 

D08  I«1.N)K) 

READ) 16,4 )COO£(K. I ) . TIME ( K , I ) , LAT) K , I ) ,LON(K, 1 ) . FLW(K. I ) , 
HSS)K, I ) ,SFW)K, I ) ,SLP(K. I ) 

FORMAT) 12 . 1X,I6,1X.F3.1.1X.F4.1,1X.I5.1X.I4,1X.I5.1X.I4) 

CONTINUE 
CONTINUE 
0077  K=1 ,NN 
0088  1=1, N)K) 

:tim=tim£)k, I ) 

CALL  MINUTE) ITIM, MIN) 

TIME)K, I )*MIN 
IF ) CODE! K . I  )  . EO.  10) 

DZ=C4-HSS)K, I ) 

SLP)K, I )  =  INT) ( 1000  /( EXP(DZ-C5) ) )* . 5  J 
ELSE 

IF)CODE)K, I ) . E0.20) 

HSS)K, I )-HSS)K, 1+1 ) 

IF)SLP)K, I ) .EO.O) 

SLP1K,I)=INT))  1 15*HSS(K, I ) )+645 . ) 

else 

IF)HSS(K,I)  EO  0 )HSS) K, I )  =  INT) ) SLPl K , I  )-645  )/ ,  1 1 5 ) 

END  IF 
ELSE 

IF)CODE)K, I )  EO  30) 

IF)SLP) K, I ) . EO  0) 

IF)CODE)K, 1-1 )  NE .20) 

SLP)K. I )*INT) )  1 15*HSS(K, I ))+645. ) 

ELSE 

IF) SLP) K, I- 1 )  GT . 800)SLP(K, I )«SLP)K, I- 1 ) 

END  IF 
END  IF 

IF)HSS)K, I ) . EO  0)HSS)K, I )*INT( (SLP)K, I )-645)/ . 115) 

ELSE 

IF )CODE )K, I) . EO .2 1 . AND .SLP)K. I ) . EO.O) 
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ELSE 

IF(CODE(K. I ) .E0.41 .ANO.CODE(K.I-I) .E0.21 . AND . SLP(K . I ) . EO . 0) 
SLP(K.I)-SLP(K.I-1) 

END  IF 
END  IF 
END  IF 
END  IF 
END  IF 
88  CONTINUE 
77  CONTINUE 

PART 

THIS  PART  SEARCHES  THE  DATA  FOR  FIX  POSITIONS  AND  TIMES.  IT  TRACKS 
BOTH  SURFACE  AND  FLIGHT  LEVEL  CENTERS  WITH  TIME.  FOR  THOSE  FIXES 
WHEN  ONLY  A  FLIGHT  LEVEL  CENTER  IS  FIXED.  THE  SURFACE  CENTER  IS 
ASSUMED  TO  BE  DIRECTLY  BELOW.  THE  SYSTEM  VELOCITY  IS  COMPUTED  FOR 
USE  IN  THE  NEXT  PART.  CENTRAL  SLP  IS  TRACKED  SO  THAT  A  CONTINUOUS 
ESTIMATE  OF  CURRENT  CENTRAL  SLP  CAN  BE  INTERPOLATED. 

VARIABLE  list********** 


TLAT  -  TEMPORARY  LATITUDE 
TLON  -  TEMPORARY  LONGITUDE 
SLAT1  -  PREVIOUS  SURFACE  FIX  LATITUDE 
SL0N1  -  PREVIOUS  SURFACE  FIX  LONGITUDE 
SMIN1  -  PREVIOUS  SURFACE  FIX  MINUTE 
SCLAT(J)  -  LATITUDE  OF  JTH  SURFACE  FIX 
SCLON(J)  -  LONGITUDE  OF  JTH  SURFACE  FIX 
SCMIN(J)  -  MINUTE  OF  JTH  SURFACE  FIX 

SCVLAT(J)  -  surface  center  VELOCITY  BASED  UPON  JTH  FIXILATITUDE 
COMPONENT ) 

SCVLON(J)  -  SURFACE  CENTER  VELOCITY  BASED  ON  JTH  FIX  (LONGITUDE 
COMPONENT ) 

SCSLP(J)  -  CENTRAL  SEA  LEVEL  PRESSURE  FOR  JTH  FIX 
SFLG  -  SURFACE  FIX  FLAG 

SVELAT  -  COMPUTED  SURFACE  CENTER  LATITUDE  VELOCITY 
SVELON  -  COMPUTED  SURFACE  CENTER  LONGITUDE  VELOCITY 
JFL6  -  flag  for  flight  LEVEL  CENTER  FIX  ONLY 
LLAT1  -  PREVIOUS  LOW  LEVEL  FIX  LATITUDE 
LL0N1  -  PREVIOUS  LOW  LEVEL  FIX  LONGITUDE 
LMINI  -  PREVIOUS  LOW  LEVEL  FIX  MINUTE 
LCLAT(J)  -  latitude  OF  JTH  LOW  LEVEL  FIX 
LCLON(J)  -  LONGITUDE  OF  JTH  LOW  LEVEL  FIX 
LCMINIJ)  -  MINUTE  OF  JTH  LOW  LEVEL  FIX 

LCVLAT(JI  -  LOW  LEVEL  CENTER  LATITUDE  VELOCITY  BASED  ON  JTH  FIX 
LCVLONIJ)  -  LOW  level  CENTER  LONGITUDE  VELOCITY  BASED  ON  JTH  FIX 
LVELAT  -  COMPUTED  LOW  LEVEL  CENTER  LATITUDE  VELOCITY 
LVELON  -  COMPUTED  LOW  LEVEL  CENTER  LONGITUDE  VELOCITY 
JJ  -  total  fixes  for  a  STORM 

SFLG*0 
DOe  K« 1 . NN 
003  1*1. NIKI 
TLAT*LAT(K. I  I 
TL0N*L0N(K. I  ) 

MIN-TIME(K. I ) 

I F ( CODE ( K . I ) . EO . 20 . OR . CODE (K. DEO. 21) 

J«  j+  1 

IF( J.EO. 1 ) 

CALL  EXCHG(TLAT.SLAT1. TLON. SL0N1. MIN. SMIN1 ) 

SCLAT( J)«SLAT1 
SCLON( J)«SL0N1 
SCMIN( J)«SMIN1 
SCSLP( J)-SLP(K.I ) 

SCVLAT( J)»0 
SCVLON( J)-0. 
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SFL6-1 

ELSE 

CALL  VEL( TLAT . TLON . MIN. SLAT  1 . SLON1 . SMIN1 , SVELAT . SVELON ) 

SCLAT( J)»SLAT1 

SCLON(vJ)-SLON1 

SCMIN( J)>SM1N1 

SCSLP(d)«SLP(K.I ) 

SCVLAT(d)«SVELAT 

SCVLON(J) 'SVELON 

SFLG-1 

END  IF 

ELSE 

IF(CODE(K.I) .E0.41 .0R.C0DE(K.I).E0.3O) 

I F ( CODE ( K . I -  1) . GT . 2 1 . OR . CODE ( K . I -  1 ) . LT . 20 ) 

d«d*1 

dFLG«1 

ELSE 

JFLG'O 

END  IF 

IF(d.EO. 1 ) 

CALL  EXCHG( TLAT. LLAT 1 , TLON. LL0N1 .MIN. LMIN1 ) 

LCLAT(d)«LLATl 

LCLON( J)<LL0N1 

LCM1N( d I'LMINI 

LCVLAT(d)»0. 

LCVL0N(d)*0. 

IF(SFLG.NE .0)GO  TO  3 

CALL  EXCHG( LLAT1 .SLAT1 .LL0N1 .SL0N1 .LMIN1 .SWIN1 ) 

SCLAT(d)-SLAT1 

SCL0N(d)-SL0N1 

SCMIN(d)«SMlN1 

SCSLP(d)»SLP(K. I ) 

SCVLAT(d)«0. 

SCVL0N(d)«0. 

SFLG«1 

ELSE 

CALL  VEL( TLAT. TLON. MIN. LLAT 1 .LL0N1 .LMIN1 .LVELAT, LVELON) 

LCLAT(d)*LLAT1 

LCL0N(d)<LL0N1 

LCMIN(d)*LMIN1 

LCVLAT(d)«LVELAT 

LCVLONi d )*LVELON 

IF(dFLG.NE , 1 )GC  TO  3 

CALL  VEL (TLAT. TLON. MIN. SLAT  1 .SL0N1 . SMIN1 . SVELAT . SVELON) 

SCLAT(d)-SLAT1 

SClON(  J  I  =  SL0M 

SCMIN(dl*SKIN1 

SCSLP( d l=SLP(K. I ) 

SCVLATl d I'SVELAT 
SCVLON( d )-SVEl0N 
SFLG«1 
END  IF 
END  IF 
END  IF 
CONTINUE 
CONTINUE 
dd'd 

IF(dd.GT  1) 

SCVLAT( 1 )»SCVLAT(2) 

SCVLON( 1 )«SCVL0N(2 ) 

LCVLAT(  1  )'LCVLAT(2) 

LCVLON( 1 )'LCVL0N(2 ) 

END  IF 
d*  1 


PART3 . . . . . . . . . 

THIS  PART  SEARCHES  THE  DATA  FOR  OBSERVATIONS  AT  A  COMMON  LEVEL 
(EITHER  700MB  OR  1500FT)  IT  WORKS  TO  DERIVE  PRESSURE/WIND  RELA- 


84 


TIONSHIPS  FOR  BOTH  THE  FLIGHT  LEVEL  AND  SURFACE.  DATA  FROM  PART  2 
IS  USED  TO  DETERMINE  THE  SURFACE  AND  FLIGHT  LEVEL  CENTER  POSITIONS 
AT  EACH  OBSERVATION  TIME.  DISTANCE  OF  THE  INDIVIDUAL  OBSERVATION 
FROM  THE  PROJECTED  CENTER  IS  USED  ALONG  WITH  THE  ANGULAR  RADIAL 
FROM  THE  CENTER  TO  THE  OBSERVATION  POSITION  TO  PLACE  THE  OBSERVA¬ 
TION  RELATIVE  TO  THE  SURFACE  AND  FLIGHT  LEVEL  CENTERS.  LOW  LEVEL 
(1500FT)  AND  700MB  OBSERVATIONS  ARE  HANDLED  SEPARATELY.  OBSERVED 
WINDS  ARE  ADJUSTED  TO  TAKE  OUT  COMPONENTS  DUE  TO  SYSTEM  MOTION. 
THESE  WINDS  ARE  THEN  TRANSPOSED  ANGULARLY  TO  A  COMMON  RADIAL  BE¬ 
TWEEN  THE  NEXT  COMMON  LEVEL  OBSERVATION  AND  THE  VORTEX  CENTER  AT 
THE  NEXT  OBSERVATION  TIME.  WIND  DIRECTION  AND.  THEREFORE.  THE  U 
(EAST-WEST)  AND  V  (NORTH-SOUTH)  COMPONENTS  ARE  ADJUSTED  FOR  THIS 
TRANSPOSITION.  AVERAGE  AND  MAXIMUM  OBSERVED  WIND  SPEEDS  ARE  THEN 
COMPUTED  BETWEEN  OBSERVATIONS  WITH  SLP  OR  700MB  HEIGHT  INFORMATION 
AS  LONG  AS  ONLY  A  MAXIMUM  OF  ONE  INTERMEDIATE  WIND  WITHOUT  PRES¬ 
SURE  INFORMATION  OCCURS  BETWEEN.  FOR  BOTH  FLIGHT  LEVEL  AND  SURFACE 
LEVEL.  A  PRESSURE  GRADIENT  IS  COMPUTED  AND  ASSIGNED  TO  THE  MID¬ 
POINT  DISTANCE  BETWEEN  OBSERVATIONS  ALONG  THE  COMMON  RADIAL. 
FACTORS  ARE  COMPUTED  WHICH  RELATE  THE  AVERAGE  WIND  TO  THE  PRES¬ 
SURE/HEIGHT  GRADIENT  AT  FLIGHT  LEVEL  AND  AT  THE  SURFACE.  THE  MAX¬ 
IMUM  WIND  BETWEEN  OBSERVATIONS  IS  ALSO  RELATED  TO  THE  SAME.  OUTPUT 
FILE  NUMBER  1  IS  CREATED  CONTAINING  THE  INPUT  DATA  WITH  TIME  CON¬ 
VERTED  TO  CONTINUOUS  MINUTE  FORMAT.  WIND  VELOCITIES  WITH  SYSTEM 
MOTION  SUBTRACTED.  AND  WITH  SLP  EXTRAPOLATED  WHEN  NECESSARY.  OUT¬ 
PUT  file  number  2  CONTAINS  DISTANCES  TO  THE  MIDPOINT  BETWEEN  OB¬ 
SERVATIONS  AND  THE  CORRESPONDING  WIND/PRESSURE  FACTORS  FOR  THAT 
POINT.  THE  RAOIUS  OF  STORM  SIZE  IS  ESTIMATED  BV  DETERMINING  A 
VALUE  FOR  the  OUTERMOST  CLOSED  ISOBAR  AND  THEN  USING  AN  ANALYTIC 
MODEL  TO  DETERMINE  ITS  RADIUS.  OUTPUT  FILE  NUMBER  3  IS  CREATED 
CONTAINING  THIS  RAOIUS  PLUS  THE  ENVIRONMENTAL  VALUES  USED  TO 
DETERMINE  THIS  RADIUS. 

VARIABLE  LIST . . . . 

LFLG  -  INDICATES  CONSECUTIVE  OBS  WITH  OR  WITHOUT  PRESSURE  DATA 

OIFTIM  -  TIME  DIFFERENCE  BETWEEN  FIX  AND  OBSERVATION 

CLAT  -  FLIGHT  LEVEL  CENTER  LATITUDE  AT  OBSERVATION  TIME 

CLON  -  flight  level  center  longitude  at  OBSERVATION  TIME 

SFCLAT  -  SURFACE  CENTER  LATITUDE  AT  OBS  TIME 

SFCLON  -  SURFACE  CENTER  LONGITUDE  AT  OBS  TIME 

DISLLA(I)  -  LATITUDE  DISTANCE  OF  OBS(K.I)  TO  FLIGHT  LEVEL  CENTER 
AT  TIME(K. I ) 

DISLLO(I)  -  LONGITUDE  DISTANCE  OF  OBS(K.I)  TO  FLIGH”  LEVEL  CENTER 
AT  TIME(K.I) 

DISSLA(I)  -  LATITUDE  DISTANCE  OF  OBS(K.I)  TO  SURFACE  CENTER  AT 
TIME(K, I ) 

DISSLO(I)  -  LONGITUDE  DISTANCE  OF  OBS(K.I)  TO  SURFACE  CENTER  AT 
TIME! K . : ) 

Cl  -  constant  FOR  DEGREE/RADIAN  RELATION 

DISTL1  -  DISTANCE  OF  OBS(K.I)  TO  FLIGHT  LEVEL  CENTER  (DEG  LAT ) 
DISTS1  -  DISTANCE  OF  OBS(K.I)  TO  SURFACE  CENTER  (DEG  LAT) 

DISTL2  -  DISTANCE  OF  OBS(K.I-I)  TO  FLIGHT  LEVEL  CENTER  (DEG  LAT) 

DISTS2  -  DISTANCE  OF  OES(K.I-I)  TO  SURFACE  CENTER  (DEG  LAT) 

DISTL3  -  DISTANCE  OF  0BS(K.I-2)  TO  FLIGHT  LEVEL  CENTER  (DEG  LAT  I 

DISTS3  -  DISTANCE  OF  0BS(K.I-2)  TO  SURFACE  CENTER  (DEG  LAT) 

ANG..1  -  RADIAL  ANGLE  BETWEEN  FLIGHT  LEVEL  CENTER  AND  OBS(K.I) 

ANGL2  -  RADIAL  ANGLE  BETWEEN  FLIGHT  LEVEL  CENTER  AND  OBS(K.I-I) 

ANGS1  -  RADIAL  ANGLE  BETWEEN  SURFACE  CENTER  AND  OBS(K.I) 

ANGS2  -  RADIAL  ANGLE  BETWEEN  SURFACE  CENTER  AND  OBS(K.I-I) 

FLANG3  -  flight  LEVEL  WIND  DIRECTION  ANGLE  FOR  0BS(K,I-2) 

FLANG2  -  FLIGHT  LEVEL  WIND  DIRECTION  ANGLE  FOR  OBS(K.I-I) 

FLANGI  -  FLIGHT  LEVEL  WIND  DIRECTION  ANGLE  FOR  OBS(K.I) 

SLANG3  -  SURFACE  WIND  DIRECTION  ANGLE  FOR  0BS(K.I-2) 

SLANG2  -  SURFACE  WIND  DIRECTION  ANGLE  FOR  OBS(K.I-l) 

SLANG)  -  SURFACE  WIND  DIRECTION  ANGLE  FOR  OBS(K.I) 

K1  -  CONSTANT  TO  SEPARATE  DIRECTION  FROM  SPEED  DATA 

K2  -  CONSTANT  TO  DETERMINE  WIND  DIRECTION  ANGLE  WITHIN  10  DEGREES 

LCKLAT  -  FLIGHT  LEVEL  CENTER  V  VELOCITY  COMPONENT  (KTS) 
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LCKLON  -  FLIGHT  LEVEL  CENTER  U  VELOCITY  COMPONENT  (KTS) 

C2  -  CONSTANT  TO  CONVERT  FROM  (OEG  LAT)/MIN  TO  KNOTS 
SCKLAT  -  SURFACE  CENTER  V  VELOCITY  COMPONENT  (KTS) 

SCKLON  -  SURFACE  CENTER  U  VELOCITY  COMPONENT  (KTS) 

FLVEL1  -  FLIGHT  LEVEL  WIND  SPEED  FOR  OBS(K.I) 

FLVEL2  -  FLIGHT  LEVEL  WIND  SPEED  FOR  OBS(K.I-I) 

FLVEL3  -  FLIGHT  LEVEL  WIND  SPEED  FOR  OBS{K.I-2) 

FUVEL1  -  FLIGHT  LEVEL  WIND  U  COMPONENT  FOR  OBS(K.I) 

FUVEL2  -  FLIGHT  LEVEL  WIND  U  COMPONENT  FOR  OBS(K.I-I) 

FUVEL3  -  FLIGHT  LEVEL  WIND  U  COMPONENT  FOR  OBS(K.I-2) 

FVVEL1  -  FLIGHT  LEVEL  WIND  V  COMPONENT  FOR  OBS(K.I) 

FVVEL2  -  FLIGHT  LEVEL  WIND  V  COMPONENT  FOR  OBS(K.I-I) 

FVVEL3  -  FLIGHT  LEVEL  WIND  V  COMPONENT  FOR  OBS(K.I-2) 

SFVEL1  -  SURFACE  WIND  SPEED  FOR  OBS(K.I) 

SFVEL2  -  SURFACE  WIND  SPEED  FOR  OBS(K.I-I) 

SFVEL3  -  SURFACE  WIND  SPEED  FOR  OBS(K.l-2) 

SUVEL1  -  SURFACE  WIND  U  COMPONENT  FOR  OBS(K,I) 

SUVEL2  -  SURFACE  WIND  U  COMPONENT  FOR  OBS(K.I-I) 

SUVEL3  -  SURFACE  WIND  U  COMPONENT  FOR  OBS(K.I-2) 

SVVEL1  -  SURFACE  WIND  V  COMPONENT  FOR  OBS(K.I) 

SVVEL2  -  SURFACE  WIND  V  COMPONENT  FOR  OBS(K.I-I) 

SVVEL3  -  SURFACE  WIND  V  COMPONENT  FOR  OBS(K.I-2) 

LANDIF  -  ANGULAR  DIFFERENCE  BETWEEN  ANGL 1  AND  ANGL2 
SANDIF  -  ANGULAR  DIFFERENCE  BETWEEN  ANGS1  AND  ANGS2 
A  -  COUNTER  OF  MISSING  FLIGHT  LEVEL  WINDS 
K3  -  CODE  FOR  MISSING  WIND 

IFLG  -  FLAG  FOR  MID  POINT  WIND  VERIFICATION 

FLUVEL  -  AVERAGE  FLIGHT  LEVEL  WIND  U  VELOCITY  COMPONENT 

FLVVEL  -  AVERAGE  FLIGHT  LEVEL  WIND  V  VELOCITY  COMPONENT 

FLMAX  -  MAX  FLIGHT  LEVEL  WIND  SPEED  BETWEEN  OBS  WITH  PRESSURE  DATA 

B  -  COUNTER  OF  MISSING  SURFACE  WINDS 

SFUVEL  -  AVERAGE  SURFACE  WIND  U  VELOCITY  COMPONENT 

SFVVEL  -  AVERAGE  SURFACE  WIND  V  VELOCITY  COMPONENT 

SFMAX  -  maximum  surface  WIND  SPEED  BETWEEN  OBS  WITH  PRESSURE  DATA 

IFLDIR  •  AVERAGE  FLIGHT  LEVEL  WIND  DIRECTION  ANGLE 

ISFOIR  -  AVERAGE  SURFACE  WIND  DIRECTION  ANGLE 

FLVEL  -  AVERAGE  FLIGHT  LEVEL  WIND  SPEED 

SFVEL  -  AVERAGE  SURFACE  WIND  SPEED 

OISTL  -  AVERAGE  DISTANCE  OF  OBS  TO  FLIGHT  LEVEL  CENTER 

OISTS  -  AVERAGE  DISTANCE  OF  OBS  TO  SURFACE  CENTER 

C3  -  CONSTANT  OF  CONVERSION  FROM  DEG  LAT  TO  NAUTICAL  MILE 

IPSCHG  -  SURFACE  PRESSURE  CHANGE  (MB)  BETWEEN  OBS  WITH  PRESSURE 

DATA 

IPLCHG  •  flight  level  PRESSURE  HEIGHT  CHANGE  (METERS)  BETWEEN  OBS 
WITH  PRESSURE  DATA 

GRADIL  -  GRADIENT  DISTANCE  AT  FLIGHT  LEVEL  BETWEEN  OBSERVATIONS 
ALONG  A  COMMON  RADIAL 

GRADIS  -  GRADIENT  DISTANCE  AT  SURFACE  BETWEEN  OBS  ALONG  A  COMMON 

radial 

PGRADL  -  PRESSURE  GRADIENT  AT  FLIGHT  LEVEL  BETWEEN  OBS  ALONG  A 
COMMON  RADIAL 

PGRADS  -  PRESSURE  GRADIENT  AT  SURFACE  BETWEEN  OBS  ALONG  A  COMMON 

radial 

IFLW3  -  temporary  FLIGHT  LEVEL  WIND 

IFLW2  -'temporary  FLIGHT  LEVEL  WIND  (OR  INTERMEDIATE  WIND) 

ISLP1  -  TEMPORARY  SLP 
ISLP2  -  temporary  SLP 
IHSS1  -  temporary  HSS 
IHSS2  -  TEMPORARY  HSS 

IPN  -  ENVIRONMENTAL  SEA  LEVEL  PRESSURE 

IPE  -  TROPICAL  CYCLONE  SEA  LEVEL  PRESSURE  LIMIT 

PWL  -  700MB  FLIGHT  LEVEL  WINO/PRESSURE  HEIGHT  GRADIENT  FACTOR  USING 
AVERAGE  WIND 

PWLM  -  700MB  FLIGHT  LEVEL  WINO/PRESSURE  HEIGHT  GRADIENT  FACTOR  USING 
MAXIMUM  WIND 

PWSS  -  SURFACE  WIND/PRESSURE  GRADIENT  FACTOR  USING  AVERAGE  WIND 
PWSSM  -  SURFACE  WINO/PRESSURE  GRADIENT  FACTOR  USING  MAXIMUM  WIND 
PWLS  -  SURFACE  WIN0/70OM6  PRESSURE  HEIGHT  GRADIENT  FACTOR  USING 


AVERAGE  WIND 

PWLSM  -  SURFACE  WIND/700MB  PRESSURE  HEIGHT  GRADIENT  FACTOR  USING 
MAXIMUM  WIND 

PWSL  -  LOW  LEVEL  WINO/SURFACE  PRESSURE  GRADIENT  FACTOR  USING 
AVERAGE  WIND 

PWSLM  -  LOW  LEVEL  WINO/SURFACE  PRESSURE  GRADIENT  FACTOR  USING 
MAXIMUM  WIND 

ISLPM  -  UPPER  LIMIT  OF  SLP  FOR  RADIAL  SIZE  COMPUTATION 
DISTM  -  DISTANCE  OF  LAST  OBS  IN  RADIAL  SIZE  COMPUTATION 

DISTN  -  OUTER  LIMIT  DISTANCE  FOR  MINIMUM  WIND  RADIAL  SIZE  COMP 

OISTO  -  INNER  LIMIT  DISTANCE  FOR  RADIAL  SIZE  COMPUTATION 

DISTP  -  INNER  LIMIT  DISTANCE  FOR  MAXIMUM  WIND  RADIAL  SIZE  COMP 

RS  -  RADIAL  SIZE 

BB  -  PARAMETER  FOR  LOG  PRESSURE  PROFILE 

N1  -  NUMBER  OF  PRESSURE  GRADIENT  COMPUTATIONS 

IPED  -  CORRECTION  TO  IPE  FOR  OIURNAL  ATMOSPHERIC  TIDE 

CC  -  TIME  DIFFERENCE  BETWEEN  OBS  AND  NEXT  SURFACE  FIX 

DD  -  TIME  DIFFERENCE  BETWEEN  SURROUNDING  SURFACE  FIXES 

RAT  -  RATIO  OF  TIME  BETWEEN  FIXES 

ISCSLP  -  INTERPOLATED  VALUE  OF  CURRENT  CENTRAL  SLP 
EE  -  DIFFERENCE  BETWEEN  ENVIRONMENTAL  AND  CENTRAL  SLP 
FF  -  DIFFERENCE  BETWEEN  OBSERVED  AND  CENTRAL  SLP 
AA  -  PARAMETER  FOR  LOG  PRESSURE  PROFILE 

GG  -  DIFFERENCE  BETWEEN  OUTER  LIMIT  ISOBAR  AND  CENTRAL  SLP 
0013  K>1 ,NN 

WRITE ( 17. lO)NUM(K) ,N(K) 

WRITE( ie, lO)NUM(K) .N(K) 

FORMAT( IX. 17, 7X. 12) 

IPE-1008 

ISLPM*2000 

DISTM«0. 

DISTN»1000. 

DIST0«0. 

DISTP-0. 

RS«0. 

SB«0. 

N1«0 

0012  I«1.N(K) 

IF(LCMIN( J) . LT .TIME(K, I ) . AND . J . LT , UJ )U«U* 1 
ISLP1«SLP(K, I ) 

IHSS1-HSS(K, I ) 

IF(CCDE(K,I ) .EC.  11  ) 

KLL-KLL+1 

1F(KLL.E0.2.AND.SLP(K.I-1).EO.O)KLL«1 
D:FTIM»LCM1N( J)-TIME(K.I ) 

CLAT«LCLAT( J)-LCVLAT( J)*DIFTIM 
CL0N*LCL0N( J)-LCVL0N(U)'DIFTIM 
DIFT1M=SCWIN(0)-TIME(K, I ) 

SFCLAT.SCLAT(0)-SCVLAT( J)*DIFTIM 
SFCLON»SCLON(d)-SCVLON( J)*DIFTIM 
DISLLA( n*LAT(K. I )-CLAT 

DISLLOl I  )•( LON(K. 1  ) -CLON ) 'COSI LAT ( K . I )/C1 ) 

DISSLA( I  )«LAT(K. I  )-SFCLAT 

DISSLO( I )*( LON( K , I )-SFCL0N)-C0S( LAT(K . I ) 'C 1) 
DISTL1*S0RT(01SLLA( I )--2*DISLL0( n-"2 ) 

DI STS  1 'SORT (OISSLA( I )--2+DISSL0( I )-*2 ) 

LCKLAT«LCVLAT(0)'C2 

LCKL0N«LCVL0N( d )• ( COS( LCLAT( d)/C1 ) )*C2 
SCKLAT*SCVLAT(d)*C2 

SCKL0N«SCVL0N(d)*(C0S(SCLAT(d)/C1 ) )*C2 
IF(KLL.EO  1 ) 

LFLG«0 

KHH*0 

ANGL1-ATAN2(DISLL0(I ).DISLLA(I ) )*C1 
IF(ANGL1  LEO  )ANGL1*360.«ANGL1 
ANGS1«ATAN2(DISSL0(I ).0ISSLA(I ))*C1 
1F(ANGS1  LE .0. )ANGS 1-360. ♦ANG51 


FLANG1»(FLW(K, I )/K1 )*K2 
SLANG1»(SFW(K. I )/K1 )*K2 

FLVEL1«FL0AT(FLW(K, I )-(FLW(K. I )/K1 )*K1 ) 

CALL  COMP ( FLVEL 1 , F LANG1 . C 1 . LCKLAT . LCKLON , FUVEL 1 . F VVEL 1 ) 
sfveli-float(sfw(k.i)-(sfw(k.i)/ki)»ki) 

CALL  COMP(SFVEL1 .SLANG1 ,C1 .SCKLAT.SCKLON.SUVEL1 ,SVVEL1 ) 
GO  TO  34 
END  IF 

IFdSLPI  .LT.800) 

IF(SLP(K.I+1 ) .LT. BOO. OR. SLP(K. 1-1). LT.800) 

FLANG1-900 
SLANG 1*900 
GO  TO  34 
END  IF 
LFLG*1 
END  I F 

IF(SLP(K,I-1).LT. 800 . AND . LFLG . EO . 0 ) 

F LANG 1*900 
SLANG 1*900 
GO  TO  34 
END  IF 

IF(LFLG.E0.2) 

DISTL3*DISTL2 
DISTS3*DISTS2 
END  IF 

DISTL2*S0RT(DISLLA(I-1)“»2+0ISLL0( 1-1 )**2) 
DISTS2*S0RT(DISSLA(I-1 )-*2+DISSL0( 1-1 )--2) 

IF(0ISTL2.NE .0. ) 

ANGL1«ATAN2(DISLL0(I).D1SLLA(1 ))-C1 
IF( ANGL1 . LE .0. ) ANGL 1 *360 . +ANGL 1 
ANGL2*ATAN2(DISLL0( 1-1 ) .DISLLA(  1-1 ) )-C1 
IF(ANGL2.LE.0. )ANGL2*360. ♦ANGL2 
END  IF 

IF(DISTS2.NE.O. ) 

ANGS1*ATAN2(0ISSL0(I ).DISSLA(I ))*C1 
IF( ANGS1 .LE .0. )ANGS1*360.*ANGS1 
ANGS2*ATAN2(DISSL0(I-1 ).DISSLA( 1-1 ) )*C1 
IF(ANGS2.LE.O. )ANGS2*360.+ANGS2 
END  IF 

IF(LFLG.E0.2) 

FLANG3*FLANG2 
SLANG3*SLANG2 
END  :* 

FLANG1*(FLW(K, I )/K1 )*K2 

SLANG1*(SFW(K, I )/K1 )-K2 

FLANG2*( FLW( K. I- 1 )/K1 )‘K2 

SLANG2»(SFW(K, 1-1 )/K1 )-K2 

FLVEL 1 =FL0AT( FLW(K. I )-( FLW(K, I )/K1 )-K1 ) 

CALL  COMP( FLVEL 1 , FLANG1 .C1 . LCKLAT . LCKLON . FUVE L 1 . FVVELI ) 
FLVEL2*FL0AT(FLW(K. I-1)-(FLW(K. 1- 1 )/K1 )*K1 ) 

CALL  COMP( FLVEL2 . FLANG2 , Cl . LCKLAT . LCKLON , FUVEL2 . F VVEL2 ) 
SFVEL1*FL0AT(SFW(K. I )-(SFW(K. I )/K1 )*K1) 

CALL  C0MP(SFVEL1 .SLANG1 .Cl .SCKLAT.SCKLON.SUVELI .SVVELI ) 
SFVEL2*FlOAT{ SFW(K. 1-1 )-(SFW(K. 1-1 )/K1 )-K1 ) 

CALL  COMP( SFVEL2 , SL ANG2 . C 1 , SCKL AT . SCKLON . SUVE L2 . S VVE L2 ) 
IF(LFLG.E0.2) 

FLVEL3*FL0AT(FLW(K.I-2)-(FLW(K,I-2)/K1 )-K1 ) 

CALL  COMP( FLVEL3 , FLANG3 . C1 . LCKLAT , LCKLON . FUVEL3 . FVVEL3 ) 
SFVEL3*FL0AT(SFW(K, I-2)-(SFW(K. 1-2 )/K1 )*K1 ) 

CALL  COMP ( SF VEL3 . SLANG3 , Cl , SCKLAT , SCKLON , 5UVEL3 . S VVEL3 ) 
END  IF 

IF(DISTL2.NE .0. ) 

LAND I F  * INT ( ANGL 1  - ANGL2 ) 

FLANG2*FLANG2*LANDIF 
IF( FLANG2 . GT . 360)FLANG2*FlANG2-360 
IF(FLANG2 . LE . 0)F LANG2*360*FLANG2 
END  IF 


SAN0IF>INT(ANSS1- ANGS2 ) 

SLANG2>SLANQ2-»SANDIF 
I F ( SLANG2 . GT . 360 ) SLANG2>SLANG2-360 
I F  ( S  L  ANG2  .  L  E  .  0 )  S  L ANG2  >3604^  SL ANG2 
END  IF 

IF(LFLG.E0.2) 

F  L  ANG3>  F  L  ANG3->'  LAND  I F 
SLANG3>SLANG3-^LAN0IF 
I F ( F  LANG3 . GT . 360 ) F  L ANG3-F  LANG3- 360 
I F ( FLANG3 . LE . 0 ) FLANG3>360*FLANG3 
I F ( SLANG3 . GT . 360 )SLANG3>SLANG3-360 
I F  ( SLANG3 .  LE  .  0 )SLANG3<>3604'SLANG3 
A«0. 

1F(FLW(K.I-2).E0.K3) 

A»A*1 . 

FUVEL3»0. 

FVVEL3*0. 

END  IF 

CALL  MI0(DISTL1 .0ISTL2.DISTL3,IFLG.FUVEL2,FVVEL2. A) 
IF(IFLG-1 )49,50.49 
IF(FLW(K. 1-1 ) .E0.K3) 

A«A*1  . 

FUVEL2«0. 

FVVEL2»0. 

END  IF 

IF(FLW(K. 1 ) .E0.K3) 

A*A* 1 . 

FUVEL1«0. 

FVVEL1*0. 

END  IF 
IF(A.E0.3.  ) 

FLUVEL»5O0. 

FLVVEL»500. 

ELSE 

FLLIVEL»(FUVEu1-*'FUVEL2*FUVEL3)/(3.0-A) 

FLVVEL«(FVVEL1*FVVEL2+FVVEL3)/<3.0-A) 

END  IF 
B«0. 

1F(SFW(K.I-2).E0.K3) 

B«B*1 . 

SUVEL3*0. 

SVVEL3«0. 

END  !'■ 

CALL  MID(DISTS1 . DI STS2 . DI STS3 . I FLG . SUVEL2 , SVVEL2 . B ) 
IF( IFLG- 1)48.51.48 
IF( SFW( K . I  -  1  )  EC  K3 ) 

B*B*1 . 

SUVEL2*0 . 

SVVEL2=0 
END  IF 

IF(SFW(K, I ) . E0.K3) 

B«B*1 . 

SUVELi'O 
SVVELI'O 
END  IF 
IF(B.E0.3  ) 

SFUVEL*50C 

SFVVEL»500 

ELSE 

SFUVEL* ( SUVEL 1+SUVEL2*SUVEL3 )/ ( 3  0-B ) 
SFVVEL*{SVVEL1*SVVEL2*SVVEL3)/(3.0-B) 

END  IF 
ELSE 
A«0  . 

IF(FLW(K. I- 1 )  E0.K3) 

A«A+1 . 

FUVEL2»0. 

FVVEL2«0. 


END  IF 

IF(FLW(K.I).E0.K3) 

A-A+1  . 

FUVEL1«0. 

FVVEL1«0. 

END  IF 
IF(A.E0.2.  ) 

FLUVEL«500. 

FLVVEL»500. 

ELSE 

FLUVEL»(FUVEL1+FUVEL2)/(2.0-A) 

FLVVEL»(FVVEL1*FVVEL2)/(2.0-A) 

END  IF 
FUVEL3-0. 

FVVEL3-0. 

B»0. 

IF(SFW(K. 1-1 ) .E0.K3) 

B-B+1  . 

SUVEL2*0. 

SVVEL2-0. 

END  IF 

IF(SFW(K. I  )  . E0.K3) 

B*B*1  . 

SUVEL 1*0 
SVVEL1-0. 

END  IF 
IF(B.E0.2.  ) 

SFUV£l*500. 

SFVVEL*500. 

ELSE 

SFUVEL- ( SUVEL 1+SUVEL2 )/( 2 .0-B ) 
SFVVEL«(SVVEL1*SVVEL2)/(2.0-B) 

END  IF 
SUVEL3«0. 

SVVEL3-0. 

END  IF 

CALL  MAXI ( FUVEL 1 . FUVEL2 . FUVEL3 . FVVEL 1 . FVVEL2 . FVVEL3 . FLMAX ) 
CALL  MAXKSUVELI , SUVEL2 . SUVEL3 . SVVEL 1 . SVVEL2 . SVVEL3 . SFMAX ) 
IFLDIR«INT( ATAN2(FLUVEL.FLVVEL)*C1-1B0. ) 

ISFDIR«INT( ATAN2(SFUVEL.SFVVEL)-C1-1B0. ) 

IF(  IFLDIR . LE .0)1FL0IR*360*IFLDIR 
IF(  ISFDIR . LE .0)ISFDIR*360*ISFDIR 
r^.VEL’^SORTlFLUVEL-  '2+FLVVEL''’2  ) 

SF VEL-SORT ( SFUVEL* •2+SFVVEL**2) 

IF(LFLG.E0.2) 

LFLG-C 

ISLP2*SLP( K. I -2  ) 

DISTL*( (DISTL3*DISTL1  )/2  .  )-C3 
DI5TS*( (DISTS3*DISTS1)/2.  )'C3 
IPSCHG*ISLP2-ISLP1 
IF( IPSCHG)47 .34.47 
GRADIL*(0ISTL3-DISTL1)*C3 
GRADIS*(DISTS3-DISTS1 )-C3 
IFIGRADIL  EO  0  ) 

PWSL«-999 

PWSLM*-999 

ELSE 

PGRADL*IPSCHG/GRADIL 
PWSL*FlVEL/PGRADL 
PWSLM'FLMAX/PGRAOL 
IF(FLVEL.GT ,500  ) 


PWSL*-999. 

PWSLM«-999. 

END  IF 
END  IF 

IF(GRADIS.E0.0  ) 
PWSS«-999. 
PWSSM--999. 


GO  TO  129 
ELSE 

PGRAOS- I PSCHG/GRAD I S 
PWSS>SFVEL/PGRAOS 
PWSSM-SFMAX/PGRAOS 
IF(SFVEL.GT.500. ) 

PWSS*-999. 

PWSSM>-g99. 

END  IF 
END  IF 
GO  TO  129 
ELSE 

IF(LFLG.EO. 1) 

LFLG-2 
GO  TO  34 
ELSE 

ISLP2-SLP1K.I-1) 
DISTL-((DISTL2*0ISTL1)/2.  )«C3 
DISTS-((0ISTS2*DISTS1)/2. )«C3 
IPSCHG-ISLP2-ISLP1 
IF( IPSCHG)44.34,44 
GRA01L*(DISTL2-DISTL1 )*C3 
GRADIS*(DISTS2-DISTS1 )*C3 
IF(GRADIL  EO.O  ) 

PWSL«-999 . 

PWSLM--999. 

ELSE 

PGRADL*IPSCHG/GRAD1L 
PWSL*FLVEL/PGRADL 
PWSLM«FLMAX/PGRADL 
IF(FLVEL.GT.500,  ) 

PWSL«-999. 

PWSLM»-999 
END  IF 
END  IF 

IF(GRAOIS.EO.O. ) 

PWSS»-999 . 

PWSSM*-999. 

GO  TO  129 
ELSE 

PGRADS'IPSCHG/GRAOIS 
PWSS-SFVEL/PGRADS 
PWSS^’-SFMA>,.'PGRADS 
IF(SFVEL  GT .500  ) 

PWSS*-999 
PWSSM»-999 
END  IF 
END  It 
GO  TO  129 
END  IF 
END  IF 
ELSE 

IF(CODE(K,I  )  E0.21  > 

BB-  5 
KlL*KLl* ’ 

IFiKLL  £0  2  AND,SLP(K.I-1  )  E0.0>KLL  =  1 

DISSLAI I )-0 

DISSlOI I  )-C 

DISlLAI I  )«0 

OISLLOI  I  I-O 

if(kll,eo  1) 

F LANG  1*900 
SLANG  1*900 
LFLG*0 
KHH*0 
GO  TO  34 
END  IF 

IF( ISLPI ,LT .800) 


tit. 


FLAN61-900 
SLANS1>900 
GO  TO  34 
END  IF 

IF(SLP(K,I-1 ) .LT.8OO.AND.LFLG.E0.O) 

FLANG1*900 
SLANG 1 >900 
GO  TO  34 
END  IF 
DISTL1-0. 

DISTS1-0. 

FLANG1»(FLW(K.I)/K1)>K2 
SLANG1»(SFW(K. I )/K1 )>K2 
IF(LFLG.E0.2) 

DISTS3-DISTS2 
SLANG3«SLANG2 
END  IF 

DISTS2>S0RT(DISSLA( I-1)-*2+DISSL0(I-1 )>*2) 

SLANG2«(SFW(K. 1-1 )/K1 )>K2 
SCKLAT>SCVLAT(d)»C2 

SCKL0N«SCVL0N(d)*(C0S(SCLAT(d)/C1 ) )»C2 
SFVEL2«FL0AT(SFW(K. 1-1 )-(SFW(K.I-1  )/K1  )*K1  ) 

CALL  COMP( SFVEL2 . SLANG2 .Cl. SCKLAT . SCKLON . SUVEL2 . SVVEL2 ) 
B*0. 

IF(LFLG.E0.2) 

SFVEL3«FL0AT(SFW(K. I-2)-(SFW(K. I-2)/K1 )-K1 ) 

CALL  COMP(  SFVEL3  .  SLANG3  .Cl.  SCKLAT  .  SCKLON .  SLIVEL3  .  SVVEL3  ) 
IF(SFW(K. 1-2) .E0.K3) 

B«B*1  . 

SUVEL3-0 

SVVEL3=0. 

END  IF 

CALL  MID(DISTS1 .0ISTS2.DISTS3.IFLG.SUVEL2.SVVEL2.B) 

IF( IFLG-1 )45.53.45 
IF(SFW(K.I-1).E0.K3) 

B>B* 1  . 

SUVEL2»0. 

SVVEL2-0. 

END  IF 
IF(B.E0.2.  ) 

SFUVEL>500. 

SFVVEL-500. 

ELSE 

SFUVEL>(SUVEL2+SUVEL3)/(2.0-B) 

SFVVEL-(SVVEL2*SVVEL3)/(2.0-B) 

SUVELI'C. 

SVVEL1=0. 

END  IF 
ELSE 

IF(SFW(K. 1-1 )  NE .K3) 

SFJVEL>SUVEL2 

SFVVEL»SVVEl2 

SUVEL1>0, 

SVVELI'O. 

SUVEL3=C 

SVVEL3=0. 

ELSE 

SFUVEL«500 

SFVVEL=50C. 

END  IF 
END  IF 

CALL  MAXKSUVELI .SUVEL2.SUVEL3.SVVEL1 .SVVEL2.SVVEL3.SFMAX) 
IFLDIR>-1 

ISFDIR-INTl AT4N2(SFUVEL.SFVVEL)*C1-180.  ) 

IF( ISFDIR. LE  0)ISFDIR*360+I5FDIR 
FLVEL*999. 

SFVEL*S0RT(SFUVEL**2*SFVVEL**2) 

1F(LFLG.E0.2) 


ISLP2-SLP(K,I-2) 

DISTS«(DISTS3/2.  )»C3 
IPSCHG-ISLP2-ISLP1 
IF(IPSCHG)140.34. 140 
GRADIS«DISTS3»C3 
PWSL*-999. 

PWSLM>-999. 

IF(GRADIS)143.34. 143 
PGRADS* I PSCHG/GRAO I S 
CALL  BEE(PGRADS.BB) 
PWSS>SFVEL/PGRADS 
PWSSM>SFMAX/PGRADS 
1F(SFVEL.GT.500. ) 

PWSS--999. 

PWSSM--999 . 

END  IF 
GO  TO  129 
ELSE 

ISLP2«SLP(K. 1-1 ) 

DISTS«(DISTS2/2 . )*C3 
IPSCHG-ISLP2-ISLP1 
IF( IPSCHG)144.34, 144 
GRADIS*DISTS2-C3 
PWSL*-999 . 

PWSLM--999 . 

1F(GRADIS)147,34, 147 
PGRADS*IPSCHG/GRADIS 
CALL  BEE(PGRAOS.BB ) 
PWSS»SFVEL/PGRADS 
PWSSM«SFMAX/PGRADS 
IF(SFVEL.GT .500  ) 

PWSS*-999. 

PWSSM«-999 . 

END  IF 
GO  TO  129 
END  IF 
ELSE 

IF(CODE(K, I ) .E0.41 ) 

KLL*KLL*1 
BB«  .  5 

IF(KLL .E0.2 . AND.SLPCK. 1-1 ) .EO.O)KLL" 
D1SLLA( I  )=0 
OlSLLOt: )=C. 

DISSLA( I )«0. 

DISSLO( I )=0. 

I'^IKLL  .EO-  1  I 
FLANG1*900 
SlANGi =300 
LFLG-0 
KHH«0 
GO  TO  34 
END  IF 

IF( ISLP1 .LT  800) 

FLANGI =900 
SLANG1=900 
L  F  LG  =  0 
GO  TO  34 
END  IF 

IF(CODE(K. 1-1 )  NE .21) 

IF( SLP(K, 1-1 ) . LT .800  AND.LFLG.EO  0) 

F LANG  1=900 
SLANG  1=900 
GO  TO  34 
END  IF 

FLANG1«(FLW(K, I )/K1 )*K2 
SLANG1=(SFW(K,I )/K1 )=K2 
IF(LFLG.E0.2) 

DISTS3=DISTS2 


.  AND.SLPCK. 1-1 ) .EO.O)KLL  = 


SLANG3«SLANG2 
SUVEL3-SUVEL2 
SVVEL3-SVVEL2 
END  IF 

DISTS2-DI3TS1 

DISTS1-0. 

SLANG2*SL ANG1 
SUVEL2>SUVEL1 
SVVEL2*SVVEL1 
B-0. 

1F(LFLG.E0.2) 

IF(SFW(K.I-2) .E0.K3) 

B“B*  1  . 

SUVEL3*0. 

SVVEL3-0. 

END  IF 

CALL  MID(DISTS1 . DI STS2 . DISTS3 . I FLG . SUVEL2 . SVVEL2 . B ) 
IF(IFLG-1 ) 205. 206, 205 
IF(SFW(K. 1-1 ) . E0.K3) 

B»B+1  . 

SUVEL2=0 

SVVEL2*0. 

END  IF 
IF(e  E0.2  ) 

SFUVEL*500. 

SFVVEL*500. 

ELSE 

SFUVEL  =  (  SUVEL2*SUVEL3  ).'(2  ^  -B  ) 

SFVVEL'(SVVEL2*SVVEL3)/(2  -B ) 

SUVEL1=0. 

SVVEL1=0. 

END  IF 
ELSE 

IF(SFW(K, 1-1 ) .NE .K3) 

SFUVEL-SUVEL2 

SFVVEL»SVVEL2 

SUV£L1*0. 

SVVEL1*0. 

SUVEL3»0. 

SVVEL3=0. 

ELSE 

SFUVEL=500. 

SFVVEL=50C . 

END  IF 
END  IF 

CALL  MAXI  (  SUVEL  1  ,  S'JVEL2  .  SUVEL3  .  SVVEL  1  ,  SVVEL2  ,  SVVEL3  ,  S'^MAX  ) 
ISFDIR=INT( ATAN2( SFUVEL.SFVVEL )-C1- 18C  ) 

IF(  ISFDIR  .  LE  •OIISFDIR'SGO-'ISFDIR 
SFVEL*SORT( SFUVEL- •2*SFVVEL- *2 ) 

IF(SFVEL.GT  500. )SFVEL*999. 

ELSE 

DISTSI'O 

IF ( SLP( K. I -2  I  LT.800  AND.LFlG.EO  0)GO  TO  3J 
END  IF 

IF  I lFlG . EO . 2 ) 

DISTL3  =  DlSTi.2 
FlANG3=FLANG2 
END  IF 
DISTL1*0. 

FLANG2-FLANG1 
IF(CODE(K. I- 1 ) .E0.21 ) 

DISTL2*S0RT(DISLLA( I-2)**2*DISLL0( I-2)**2 ) 

IFLW3-FLW(K. 1-3) 

IFLW2*FLW(K,I-2) 

ELSE 

DISTL2-S0RT(DISLLA( 1-1 )’-2*DISLL0( 1-1 )--2 ) 

IFLW3*FLW(K. I -2 ) 

IFLW2»FLW(K.I-1 ) 


END  IF 

IF(LFLG.E0.2) 

FUVEL3-FUVEL2 
FVVEL3-FVVEL2 
END  IF 

FUVEL2-FUVEL1 
FVVEL2«FVVEL1 
A«0 . 

IF(LFLG.E0.2) 

IF(IFLW3.E0.K3) 

A-A+1 . 

FUVEL3-0. 

FVVEL3-0. 

END  IF 

CALL  MID(DISTL1 .DISTL2.DISTL3. IFLG.FUVEL2.FVVEL2, A) 

IF(IFLG-1 )251 ,252,251 

IF(IFLW2.E0.K3) 

A*A+1 . 

FUVEL2=0. 

FVVEL2=0. 

END  IF 
IF(A.E0.2.  ) 

FLUVEL«500. 

FLVVEL«500. 

ELSE 

FLUVEL-( FUVEL2+FUVEL3 )/( 2 -O-A ) 

FLVVEL*(FVVEL2+FVVEL3)/(2 .0-A) 

FUVEL1*0. 

FVVEL1=0. 

END  IF 
ELSE 

IF(IFLW2.NE.K3) 

FLUV£L»FUVEL2 

FLVVEL«FVVEL2 

FUVEL1*0. 

FVVEL1-0. 

FUVEL3*0. 

FVVEL3»0. 

ELSE 

FLUVEL«500. 

FLVVEL=500. 

END  IF 
END  IF 

CALL  MAXI ( FUVEL1 , FUVEL2 , FUVEL3 , FVVEL 1 , F VVEL2 , F VVEL3 . FLMAX ) 
IFLDIR-INT{ ATAN2(FLUVEL.FLVVEL)-C1- 180. ) 

IFICODEIK, I- 1 ) . EO. 21  ) 

SFVEL=999 
ISFDIR'-  1 
END  IF 

IF( IFlDIR. LE .0)IFLDIR=360+IFLDIR 
FLVEL=S0RT(FLUVEL**2+FLVVEL**2) 

IF(LFLG  E0.2) 

lflg=o 

DISTL*(DISTl3/2 ,  )*C3 
IF( CODE(K. I-  1  )  . E0.21  ) 

ISLP2*SLP( K, 1-3  ) 

ELSE 

ISLP2»SLP(K. 1-2 ) 

END  IF 

IPSCHG*ISLP2-ISLP1 

IF(IPSCHG)289,34,289 

GRADIL-DISTL3*C3 

IF(GRADIL)288,34.288 

PiSRADL'IPSCHG/GRADIL 

PWSL*FLVEL/PGRADL 

PWSLM«FLMAX/PGRADL 

IF(FLVEL .GT .500. ) 

PWSL--999. 


PWSLM--999. 

END  IF 

IF{C0DE(K,I-1).E0.21) 

PWSS>-999. 

PWSSM--999. 

P6RADS>-9.9999 

ELSE 

PGRADS'PGRADL 
CALL  BEE(PGRADS.BB) 

IF(SFVEL.E0.999. ) 

PWSS--999. 

PWSSM--999. 

ELSE 

PWSS«SFVEL/PGRAOS 

PWSSM«SFMAX/PGRADS 

END  IF 

END  IF 

GO  TO  129 

ELSE 

DISTL»(DISTL2/2. )»C3 
IF(COOE(K. 1-1 ) .E0.21) 
lSLP2»SLP(K.I-2) 

ELSE 

ISLP2«SLP(K. 1-1 ) 

END  IF 

IPSCHG«1SLP2-ISLP1 

1F(IPSCHG)244.34.244 

GRAD1L«DISTL2»C3 

IF(GRADIL)240.34.240 

PGRA0L*IPSCHG/GRAD1L 

PWSL=FLVEL/PGRADL 

PWSLM«FLMAX/PGRAOL 

IF(FLVEL.GT.500. ) 

PWSL»-999. 

PWSLM--999. 

END  IF 

IF(COOE(K. 1-1 ) .E0.21 ) 

PWSS«-999. 

PWSSM«-999. 

PGRA0S»-9.9999 

ELSE 

PGRA0S«PGRADL 
CALL  BEE(PGRADS,BE ) 

IF(SFVEL.E0.999.  ) 

PWSS»-999. 

PWSSM«-999 
GO  TO  129 
ELSE 

PWSS*SFVEL/PGRADS 

PWSSM-SFMAX/PGRADS 

END  IF 

END  IF 

END  IF 

GO  TO  129 

END  IF 

END  IF 

END  IF 

GO  TO  1  10 

N1«N1*1 

PWL*-999 . 

PWLM>-999 . 

PWLS*-999. 

PWLSM>-999. 

WRITE! 18. 123 )DI STL , OISTS . PWL . PWSS , PWLM , PWSSM , PWLS . PWLSM , 
PWSL.PWSLM 

FORMAT! 2( 1X.F4 .0) .8! 1X.F6.0) ) 

WRITE! 17,5)C0DE!K.I ).TIME!K.I ).LAT!K.I ).l0N!K,I ).FLANG1, 
FLVELI .HSS!K, I ) .SLANGI .SFVEL1 .SLP!K, I ) 
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FORMAT ( 1X.:2.1X.I6.1X.F4.1,1X.F5.1.2(1X,I3.'/'.F4.0.1X.I4)) 
60  TO  196 

IF(C0DE(K.I).E0.1O) 

KHH«KHH* 1 

I F ( KHH . EO . 2 . AND .HSS(K.I-1).E0.0 )KHH» 1 
DIFTIM-LCMIN(d)-TIME(K.I ) 

CLAT-LCLAT(d)-LCVLAT( J)*OIFTrM 
CLON«LCLON(d)-LCVLON(d)*DIFTIM 
DIFTIM»SCMIN(>J)-TIME(K.l  ) 

SFCLAT-SCLAT ( J ) -SCVLAT ( d ) *01 FTIM 
SFCLON-SCLON(d)-SCVLON{d)*DIFTIM 
DISLLA( I  )«LAT(K, 1 )-CLAT 

DISLLOd )»(LON(K, I )-CLON)»COS(LAT(K. I )/C1) 

DISSLA( I )»LAT(K. I )-SFCLAT 

DISSLOd  )»(LON((<.I  )-SFCLON)*COS(LAT(K.l  )/C1) 
DISTL1«S0RT(DISLLA( I )**2+DISLL0(  I )*»2) 

DlST51«S0RT(DISSLAd  )**2+DlSSL0d)*-2) 

LCKLAT«LCVLAT(d)*C2 

LCKL0N«LCVL0N(d)*(C0S(LCLAT(d)/C1 ) )»C2 
SCKLAT«SCVLAT ( d ) *02 

SCKL0N«SCVL0N(d)*(C0S(SCLAT(d)/C1 ) )*C2 
IF(KHH. EO. 1 ) 

KLL«0 

LFLG*0 

ANGL1-ATAN2(DISLL0d  ).DISLLAd  ))-C1 
1F(ANGL1  .LE.O.  )  ANGL  1  «360 . -^ANGL 1 
ANGS1«ATAN2(DISSL0( I ).DISSLA( I ) )*C1 
IF( ANGSI .LE.O. I ANGS 1 *360 . ♦ANGS 1 
FlANG1*(FLW(K. I)/K1 )-K2 
SLANG1«(SFW(K, I )/K1 )*K2 

FLVEU1«FL0AT(FLW(K.I )-(FLW(K. I )/K1 j-KD 

CALL  C0MP(FLVEL1 . F LANG  1 . C 1 . LCKLAT . LCKLON . FUVEL I , FVVEL1 ) 

SFVEL1«FLQAT(SFW(K.I )-(SFW(K.I )/K1 )»K1 ) 

CALL  C0MP(SFVEL1 ,SLANG1 .Cl . SCKLAT . SCKLON , SUVEL 1 .SVVEL1 ) 

GO  TO  134 
END  IF 

IFdHSSI  .EO.O) 

IF(HSS(K. 1*1 ) .EO.O.OR.HSS(K. 1-1 ) .EO.O) 

F LANG 1*900 
SLANG  1*900 
GO  TO  134 
END  IF 
LFLG=1 
END  IF 

IF(HSS(K.I-1 ) . EO.O. AND. LFLG. EO.O) 

FLANG1*900 
SLANG  1*900 
GO  TO  134 
END  IF 

DISTL1*S0RT{DISLLAd )*-2+DISLL0d )--2) 

DISTS1*S0RT(DISSLAd  )-*2+0ISSL0(  n*-2) 

IF(LFLG.E0.2) 

DISTL3*DISTL2 
DISTS3*DISTS2 
END  IF 

DISTL2*S0RT(DISLLAd-1  )  -  •  2*01  SLLDI  I  -  1  »"2  ) 
DISTS2*S0RT(DISSLA( 1-1 )-'2+DISSL0d- 1 )-*2 ) 

IF(DISTl2 .NE .0. ) 

ANGL1*ATAN2(DISLL0( I ) .DISLLAd ) )*C1 
IF( ANGL1 .LE.O. ) ANGL 1 *360 . ♦ANGL 1 
ANGL2*ATAN2(DISLL0d-1  ),DISLLAd-1  ))-C1 
IF(ANGL2.LE.O. )ANGL2*360. «ANGL2 
END  IF 

IF(DISTS2.NE.O. ) 

ANGS1*ATAN2(0ISSL0< I ).DISSLAd ) )*C1 
IF(ANGS1.LE  0. ) ANGS 1*360. ♦ANGSI 
AN6S2*ATAN2(DISSL0d-1).DISSLAd-1))*C1 
IF(ANGS2.LE.O. )ANGS2*360 . ♦ANGS2 


W>V>  f 


97 


f  mT  .•  •  .* 


END  IF 

IF(LFLG.E0.2) 

FLANG3>FLANG2 
SLANG3>SLANG2 
END  IF 

FLANG1*(FLW(K,I )/K1 )»K2 
SLANG1«(SFW(K, I )/K1)»K2 
FLANG2«(FLW(K, I-1)/Ki)»K2 
SLANG2«(SFW(K.I-1 )/K1)»K2 
LCKLAT-LCVLAT ( d ) *02 

LCKLON»LCVLON(d)«(COS(LCLAT(d)/C1 ) )»C2 
SCKLAT-SCVLAT ( d ) •C2 

SCKL0N»SCVL0N(d)»(C0S(SCLAT(d)/C1 ) )»C2 
FLVEL1»FLOAT(FLW(K.I )-{FLW(K,I)/K1 )»K1  ) 

CALL  COMP( FLVEL 1 . FLANG 1 . C1 . LCKLAT . LCKLON . FUVEL 1 . FVVEL 1  ) 
FLVEL2*FL0AT{ FLW(K, 1-1 )-(FLW(K. 1-1 )/K1 )*K1) 

CALL  COMP( FLVEL2 . FLANG2 , Cl . LCKLAT . LCKLON . FUVEL2 . F VVEL2 ) 
SFVEL1-FL0AT(SFW(K.I )-(SFW(K.I )/K1 )»K1 ) 

CALL  C0MP(SFVEL1 .SLAN61.C1 .SCKLAT.SCKL0N.SUVEL1 .SVVEL1  ) 
SFVEL2«FL0AT(SFW(K.I-1 )-(SFW(K.I-1 )/K1 )»K1  ) 

CALL  COMP( SFVEL2 . SLANG2 , Cl . SCKLAT . SCKLON . SUVEL2 . SVVEL2 ) 
IF(LFLG.E0.2) 

FLVEL3*FL0AT( FLW(K. 1-2)-(FLW(K. 1-2 )/K1  )»K1  ) 

CALL  C0MP(FLVEL3.FLANG3. Cl .LCKLAT , LCKLON , FUVEL3 , FVVEL3 ) 
SFVEL3-FL0AT(SFW(K. I-2)-(SFW(K. I-2)/K1 )»K1 ) 

CALL  COMP( SF VEL3 , SLANG3 .Cl. SCKLAT . SCKLON . SUVEL3 . SVVEL3 ) 
END  IF 

IF(D1STL2 .NE .0. ) 

LANDIF-INT( ANGL1-ANGL2) 

FLANG2<FLANG2-»LAN0IF 
IF( FLANG2 .GT . 360 ) FLANG2«FLANG2-360 
I F ( F  L ANG2 . L  E . 0 ) F  L ANG2 • 360+F  L ANG2 
END  IF 

IF(OISTS2.NE.O. ) 

SANOI F  « INT ( ANGS 1 - ANGS2 ) 

SLANG2>SLANG2«SANDIF 
I F ( SLANG2 . GT . 360)SLANG2«SLANG2-360 
I F  (  SLANG2  .LEO)  SLANG2>360-»SLANG2 
END  IF 

IF(LFLG.E0.2) 

FLANG3»FLANG3*LANDIF 
SLANG3-SLANG3-^LANDI  F 
I F { FLANG3 . GT . 360 ) r LANG3=FLANG3-36C 
I F ( FLANG3  LEO ) FLANG3»360*FLANG3 
IF ( SLANG3 . GT . 360 ) SLANG3»SLANG3-360 
I F ( SL ANG3  LEO ) SL ANG3*360*SL ANG3 
A-0. 

IF( FLW(K. 1-2 ) . E0.K3) 

A«A*  1  . 

FUVEL3-0. 

FVVEL3-0. 

END  IF 

CALL  MID (DI STL  1 . DI STL2 . DI STl3 . 1 FLG . FUVE L2 . F VVEL2 . A ) 

IF( IFLG- 1 )311 . 350.311 
IFfFLWlK. I  -  1  I . EO  K3 I 
A*A*  1  . 

FUVEL2=0. 

FVVEL2»0 
END  IF 

IF(FLW(K.I )  EO  K3) 

A»A*1 . 

FUVEL 1*0 
FVVEL1*0 
END  IF 
IF(A.E0.3.  ) 

FLUVEL-500 

FLVVEL*500. 

ELSE 


^ . 


y:-. 


v:* 


m 

^>3 


FLUVEL»(FUVEL1+FUVEL2+FUVEL3)/(3.0-A) 

FLVVEL»(FVVEL1+FVVEL2+FVVEL3)/(3.0-A) 

END  IF 
B-0. 

IF(SFW(K.I-2) .E0.K3) 

B-B+1 . 

SUVEL3-0. 

SVVEL3-0. 

END  IF 

CALL  MID(DISTS1 ,DISTS2.DISTS3.IFLG.SUVEL2.SVVEL2.B) 
IF( IFLG-1 )302.351 ,302 
IF(SFW(K. 1-1 ) .E0.K3) 

B-B+1 . 

SUVEL2>0. 

SVVEL2»0. 

END  IF 

IF(SFW(K. I ) .E0.K3) 

B-B+1 . 

SUVEL1«0. 

SVVEL1-0. 

END  IF 
IF(6.E0.3. ) 

SFUVEL«500. 

S^VVEL«500. 

ELSE 

SFUVEL«(SUVEL1+SUVEL2+SUVEL3)/(3.0-B) 

SFVVEL«(SVVEL1*SVVEL2*SVV£L3)/(3.0-B) 

END  IF 
ELSE 
A«0 . 

IF( FLW(K. :-1 ) . E0.K3) 

A«A* 1 . 

FUVEL2-0. 

FVVEL2-0. 

END  IF 

IF(FLW(K.n  .E0.K3) 

A«A*  1  . 

FUVEL1*0. 

FVVEL1«0. 

END  IF 
IF(A.E0.2. ) 

FLUVEL«500. 

flvvel»eoo. 

ELSE 

FLUVEL«(FUVEL1*FUVEL2)/(2.0-A) 

FLVVEL'( fvVEL1*FVVEL2)/(2  O-A) 

FUVEL3*0. 

FVVEL3»0 
END  IF 
B»0. 

IF(SFW(K.I-1).E0.K3) 

B«B*  1 
SUVEL2»C 
SVVEL2*C 
END  IF 

IF( SFW(K, I  )  . E0  .K3  ) 

B*B*  1 
SUVEL 1*0 
SVVELI'O 
END  IF 
IF(B.E0.3. ) 

SFUVEL-500. 

SFVVEL-500. 

ELSE 

SFUVEL-(SUVEL1*SUVEL2)/(2.0-B) 

SFVVEL-(SVVEL1+SVVEL2)/(2.0-B) 

SUVEL3>0. 

SVVEL3-0. 


r. 

I 
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END  IF 
END  IF 

CALL  MAXI ( FUVEL 1 . FUVEL2 . FUVEL3 . F VVEL 1 , F VVEL2 . F VVEL3 . FLMAX ) 
CALL  MAXKSUVELI ,SUVEL2.SUVEL3,SVVEL1 ,SVVEL2.SVVEL3.SFMAX) 
IFLDIR»INT(ATAN2(FLUVEL.FLVVEL)«C1-180. ) 
ISFDIR»INT(ATAN2(SFUVEL.SFVVEL)»C1-180. ) 

IFCIFLDIR.LE .0)IFLDIR«36CHIFLDIR 
IF(ISFDIR.LE.0)ISF0IR>36CHISF0IR 
FLVEL-S0RT(FLUVEL**2+FLVVEL»»2) 
SFVEL-S0RT(SFUVEL»»2*SFVVEL»*2) 

IF(LFLG.E0.2) 

LFLG-0 

IHSS2>HSS(K. 1-2) 

ISLP2-SLP(K. 1-2) 

D1STL«( (DISTL3+D1STL1 )/2.  )*C3 
D1STS«( (DISTS3+D1STS1 )/2.  )-C3 
IPLCHG»1HSS2-IHSS1 
IPSCHG»ISLP2-ISLP1 
IF(IPLCHG)399. 134.399 
399  GRA0IL«(0ISTL3-0ISTL1 )«C3 

GRADIS»(0ISTS3-D1STS1 )*C3 
IF(GRA01L.EO.O. ) 

PWL«-999 . 

PWLS*-999 

PWLM«-g99. 

PWLSM>-999 . 

ELSE 

PGR ADL  *  I PLCHG/GR AD  I L 
PWL-FLVEL/PGRADL 
PWLM«FLMAX/PGRA0L 
IF(FLVEL.GT.500. ) 

PWL--999. 

PWLM>-999. 

END  IF 

PWLS-SFVEL/PGRADL 
PWLSM-SFMAX/PGRAOL 
IF(SFVEL.GT.500. ) 

PWLS»-999. 

PWLSM--999 . 

END  IF 
END  IF 

IFCIPSCHG.EO.O.OR.GRADIS.EO.O. ) 

PWSS=-999 . 

PWSSM«-999. 

PGRADS«-9.9999 
GO  TO  229 
ELSE 

PGRADS«IPSCHG/GRADIS 
PWSS-SFVEL/PGRADS 
PWSSM«SFMAX/PGRADS 
IFISFVEL.GT.SOO. ) 

PWSS»-999. 

PWSSM«-999 . 

END  IF 
END  IF 
GO  TO  229 

else 

IFILFLG.EO  1) 

LFLG*2 
GO  TO  134 
ELSE 

IHSS2«HSS(K.I-1  ) 

ISLP2*SLP(K,I-1) 

DISTL«((DISTL2*DISTL1)/2.  )-C3 
DISTS*( (OISTS2*DISTS1 )/2 . )*C3 
IPLCHG-IHSS2-IHSS1 
IPSCHG-ISLP2-ISLP1 
IF( IPLCHG)344. 134.344 


a 


'*>  V  V  V  •  •  ‘  8*^  8*^  m"- 


GRADIL«(DISTL2-DISTL1 )»C3 
GRADIS«(DlSTS2-0ISTSi)»C3 
IF(GRAOIL.EO.O. ) 

PWL--999. 

PWLS>-999. 

PWLM--999. 

PWLSM>-999. 

ELSE 

PGRADL  « I PLCHG/GR AO I L 
PWL'FLVEL/PGRADL 
PWLM>FLMAX/PGRADL 
IF(FLVEL.GT.500. ) 

PWL--999. 

PWLM«-999 . 

END  IF 

PWLS-SFVEL/PGRADL 
PWLSM>SFMAX/PGRAOL 
1F(SFVEL.GT.500. ) 

PWLS--999. 

PWLSM>-999. 

END  IF 
END  IF 

IFdPSCHG.EO.O.OR.GRADIS.EO.O.  ) 
PWSS*-999 . 

PWSSM*-999. 

PGRA0S>-9.9999 
GO  TO  229 
ELSE 

PGHADS*IPSCHG/GRADIS 
PWSS*SFVEL/PGRADS 
PWSSMcSFMAX/PGRADS 
IF(SFVEL.GT.500. ) 

PWSS--999. 

PWSSM>-999. 

END  IF 
END  IF 
END  IF 
60  TO  229 
END  IF 
ELSE 

IF(CODE(K. I ) . E0.20) 

BB*  .5 
KHH=KHH+ 1 

IF<KHH. EO. 2 . AND .HSS(K. 1-1 ) . E0.0)KHH-1 
DISLLA( I )«0. 

DISLLO( I  )*0. 

DISSLAC I )*0. 

DISSLO( I )»0. 

IF(KHH.E0. 1 ) 

F LANG  1 *900 
SLANG1>900 
LFLG»0 
KLL*C 
GO  TO  134 
END  IF 

IF( ISLP1 . LT .800) 

FLANG1*900 
SLANG  1*900 
GO  TO  134 
END  IF 

IF(SLP(K. 1-1 )  LT .800  AND. LFLG.EO.O) 

F LANG  1*900 
SLANG  1*900 
GO  TO  134 
END  IF 
DISTL1*0. 

DISTS1-0. 

FLAN61*(FLW(K. I )/K1 )*K2 


SLANG1»(SFW(K,1)/K1)»K2 

IF(LFLG.E0.2) 

DISTS3-DISTS2 
SLANG3«SLANG2 
END  IF 

DISTS2«S0RT(DISSLA(I-1)»»2+DISSL0(I-1 )»»2) 

SLANG2»(SFW(K.1-1)/K1)»K2 

SCKLAT-SCVLAT ( J ) *02 

SCKL0N»SCVL0N( d ) • ( C0S( SCLAT( d )/C 1) ) *02 
SFVEL2»FL0AT(SFW(K.1-1)-(SFW(K.I-1 )/K1 )»K1 ) 

CALL  COMP ( SF VEL2 . SLANG2 .Cl. SCKLAT . SCKLON . SUVEL2 . S VVEL2 ) 
B-0. 

IF(LFLG.E0.2) 

SFVEL3-FL0AT(SFW(K. l-2)-(SFW(K.l-2)/K1 )»K1 ) 

CALL  COMP( SFVEL3 . SLANG3 .Cl. SCKLAT . SCKLON . SUVEL3 . SVVEL3 ) 
IF(SFW(K,1-2).E0.K3) 

B*B+1 . 

SUVEL3-0. 

SVVEL3-0. 

END  IF 

CALL  MID(0ISTS1 .DISTS2.DISTS3. IFLG.SUVEL2.SVVEL2.B) 

IF( IFLG-1 )401 .453.401 
IF(SFW(K. 1-1 ) .E0.K3) 

B«B+1  . 

SUVEL2-0. 

SVVEL2-0. 

END  IF 
IF(B.E0.2.  ) 

SFUVEL*500. 

SFVVEL'500. 

ELSE 

SFUVEL=(SUVEL2+SUVEL3)/(2.0-B) 

SFVVEL»(SVVEL2+SVVEL3)/(2.0-B) 

SUVEL1»0. 

SVVEL1«0. 

END  IF 
ELSE 

IF(SFW(K. 1-1 ) .NE .K3) 

SFUVEL-SUVEL2 

SFVVEL»SVVEL2 

SUVEL1«0. 

SVVEL1«0. 

SdVEL3=0. 

SVVEL3«0. 

ELSE 

SFUVEL*500. 

SFVVEL=500. 

END  IF 
END  IF 

CALL  MAXKSUVELI  .  SUVEL2  .  SUVEL3 .  SVVEL  1  .  SVVEL2  .  SVVEL3  .  SFMAX  ) 
IFLDIR--1 

ISFDIR«INT( ATAN2{SFUVEL.SFVVEL)*C1-1B0. ) 

IF( ISFDIR . LE .0)ISFDIR»360*ISFDIR 
FLVEL=999. 

SFVEL-SORT( SFUVEL'-2*SFVVEL»*2 ) 

:f(lflg,eo.2) 

ISLP2*SLP(K. 1-2 ) 

DlSTS-(DISTS3/2. )*C3 
IPSCHG»ISLP2-ISLP1 
IF( IPSCHG)440, 134,440 
GRADIS-DISTS3*C3 
PWL>-999. 

PWLS»-999. 

PWLM>-999 . 

PWLSM«-999. 

IF(GRADIS)443. 134,443 
PGRADS- I PSCHG/QR AD I S 
CALL  BEE(PGRADS.BB) 


PWSS-SFVEL/P6RA0S 
PWSSM*SFMAX/PGRADS 
IF(SFVEL.6T.500. ) 

PWSS«-999. 

PWSSM--999. 

END  IF 
60  TO  229 
ELSE 

ISLP2«SLP(K. 1-1 ) 

DISTS«(DISTS2/2. )»C3 
IPSCHG«1SLP2-ISLP1 
IF(IPSCHG)444, 134,444 
GRADIS-0ISTS2*C3 
PWL--999. 

PWLS«-999. 

PW,.M>-999 . 

PWLSM«-999. 

IF(GRA0IS)447. 134,447 
PGRADS>IPSCHG/GRAOIS 
CALL  BEE(PGRADS,BB) 

PWSS«SFVEL/PGRADS 
PWSSM>SFMAX/PQRAOS 
IF(SFVEL.6T.500. ) 

PWSS--999 . 
pwssM«-ggg. 

END  IF 
GO  TO  259 
END  IF 
ELSE 
B6<‘  .5 
KHH=KHH* 1 

I F ( KMH . EO . 2 . AND . HSS ( K , I -  1 ) . £0 . 0 )KHH. 1 
DISLLA(I)-0. 

DISLL0(I)«O. 

DISSLA(I)»0. 

DISSL0(n»O. 

IF(KHH.EO. 1 ) 

FLANG1«900 
SLANG  1 >900 
LFLG-0 
KLL*0 
GO  TO  134 
END  IF 

IF( IHSS1 . EO.O. AND. CODE(K. 1-1 ) .£0.20) IHSS1»HSS(K, 1-1) 
IF(IHSS1 .EO.O) 

F LANG  1=900 
SL4NG1=900 

lflg=o 

GO  TO  134 
END  IF 

IF(C0DE(K,I-1).NE.20) 

IF(HSS(K, 1-1 ) . EO.O. AND. LFLG. EO.O) 

F LANG  1=900 
SLANG  1=900 
GO  TO  134 
END  IF 

FLANG1=( FLW(K. I )/K1 )*K2 
SLANG1=(SFW(K, I )/K1 )*K2 
IF(LFLG.E0.2) 

DISTS3-DISTS2 
SLANG3=SLANG2 
SUVEL3=SUVEL2 
SVVEL3=SVVEL2 
END  IF 

DISTS2-DISTS1 

DISTS1-0. 

SLANG2=SLANQ1 

SUVEL2-SUVEL1 


V  V 


IF(LFLG.E0.2) 

1F(IFLW3.E0.K3) 

A»A*1 . 

FUVEL3-0. 

FVVEL3-0. 

END  IF 

CALL  MID(DISTL1 .D1STL2.D1STL3 . IFLG, FUVEL2 . FVVEL2 . A ) 

1F( IFLG-1 )502.552.502 
IF(IFLW2.E0.K3) 

A»A+1 . 

FUVEL2-0. 

FVVEL2-0. 

END  IF 
IF(A.E0.2.  ) 

FLUVEL-500. 

FLVVEL-5CX3. 

ELSE 

FLUVEL«(FUVEL2+FUVEL3)/(2.0-A) 

FLVVEL»(FVVEL2+FVVEL3)/(2.0-A) 

FUVELTO. 

FVVEL1-0. 

END  IF 
ELSE 

IF( IFLW2 .NE .K3» 

FLUVEL«FUVEL2 

FLVVEL-FVVEL2 

FUVEL1«0. 

FVVEL1«0. 

FUVEL3»0. 

FVVEL3»0. 

ELSE 

FLUVEL*500. 

FLWEL'SOO. 

END  IF 
END  IF 

CALL  MAXKFUVELI , FUVEL2 , FUVEL3 , FVVEL 1 . FVVEL2 . F VVEL3 . FLMAX ) 
IFLDIR»INT( ATAN2<FLUVEL,FLVVEL)*C1-180. ) 
IF(C0DE(K,I-1).EO.2O) 

SFVEL«999. 

ISFDIR«-1 
END  IF 

IF( IFLDIR.LE .0)IFLDIR*360+IFLDIR 
FLVEL*SORT( FLUVEl**2+FLVVEL‘-2) 

IF(LFLG.E0.2) 

LFLG-0 

DISTL-<0ISTL3/2 . )*C3 
DISTS-DISTL 
IF(CODE(K, 1-1  )  . EO. 20) 

IHSS2«HSS(K, 1-3 ) 

ELSE 

IHSS2«HSS(K. 1-2) 

ISLP2-SLP(K. 1-2) 

END  IF 

IPLCHG=IHSS2-IHSS1 
IF( IPLCHG)589. 134.589 
GRADIL=DISTL3'C3 
IF(GRA0IL)588. 134.588 
PGRADL«IPLCHG/GRADIL 
PWL-FLVEL/PGRADL 
PWLM'FLMAX/PGRADL 
IF(FLVEL.GT.500. ) 

PWL--999 . 

PWLM>-999. 

END  IF 

PWLS«SFVEL/PGRADL 
PWLSM»SFMAX/PGRADL 
IFCSFVEL.GT.SOO. ) 

PWLS--999. 


PWLSM>-999. 

END  IE 

IPSCH6-ISLP2-1SLP1 

IF(IPSCH6.E0.O.0R.C0DE(K.I-1).E0.2O) 

PWSS--999. 

PWSSM«-999. 

PGRA0S*-9.9999 

ELSE 

P6RA0S> I PSCHG/QR AO I L 
CALL  BEE(PGRAOS.BB) 

PWSS-SFVEL/PGRADS 
PWSSM«SFMAX/PGRADS 
1 F ( SFVEL-999 )229 . 575 . 229 
PWSS*-999. 

PWSSM--999. 

END  IF 
GO  TO  229 
ELSE 

DISTL-{DISTL2/2. )»C3 
DISTS-DISTL 
IF(CODE(K. 1-1 ) .E0.20) 

IHSS2-HSS(K. 1-2) 

ELSE 

IHSS2«HSS(K. 1-1 ) 

ISLP2*SLP(K,1-1) 

END  IF 

IPLCHG»IHSS2-IHSS1 
IF( IPLCHG)554 , 134 . 554 
GRADIL-DISTl2"C3 
IF(GRA0IL)582. 134,582 
PGRADL* I PLCHG/GRADI L 
PWL«FLVEL/PGRADL 
PWLMiFLMAX/PGRAOL 
IF(FLVEL.GT.500. ) 

PWL«-999. 

PWLM»-999. 

END  IF 

PWLSM-SFMAX/PGRADL 
PWLS«SFVEL/PGRADL 
IF(SFVEL.GT,500. ) 

PWLS*-999. 

PWLSM«-999 . 

END  IF 

IPSCHG»ISLP2-ISLP1 

IF( IPSCHG.EO.O.OR.CODE(K.I-1 ) .E0.20) 

PWSS*-999 

PWSSM«-999 

PGRADS'-9.9999 

ELSE 

PGRADS' I PSCHG/GRADI L 
CALL  BEE(PGRAOS.BB) 

PWSS*SFVEL/PGRADS 
PWSSM=SFMAX/PGRADS 
IF(SFVEL-999  1229,576.229 
PWSS=-999 . 

PWSSM=-999 
END  IF 
GO  TO  229 
END  IF 
END  IF 
END  IF 
GO  TO  12 
N1*N1*1 
PWSL«-999 . 

PWSLM--999. 

WRITE( 18. 123)0ISTL.DISTS.PWL,PWSS,PWLM.PWSSM,PWLS,PWLSM. 
PWSL.PWSLM 

WRITE( 17.5)C0DE(K,I).TIME(K,I).LAT(K.I),L0N(K.I),FLANG1, 


FLVEL1 .HSS(K.I ).SLAN61 .SFVEL1.SLP(K. I ) 

96  IF(ISLP1-800)12. 199, 199 

99  IF(SLANG1 .NE . 900. OR . FLANG1 .NE . 900) 

IFdSLPI  .GT.IPE.AND.D1STS1  .6T.DISTP) 

IF(FLVEL1 .GE.35. . OR . SFVEL 1 .GE . 25 . ) 

IPE-ISLP1 

DISTM-DISTS1 

DISTP«DISTS1 

ISLPM«2000 

END  IF 

ELSE 

IF(FLVEL1 .GE.35. .OR. SFVEL1 .GE .25 . ) 

IF(DISTS1 .GT.DISTP)D1STP-DISTS1 
ELSE 

IF(DISTS1 .LT.DISTN.AND.DISTS1 .GT .DISTP) 

IF(FLVEL1 .LT.25.  . AND . SFVEL 1 . LT . 20 . ) 

IPE*ISLP1 
D1STM-D1STS1 
DISTN-DISTS1 
ISLPM-2000 
END  IF 
END  IF 
END  IF 
END  IF 
END  IF 

IF(BB.GT.O. .AND.DISTS1 .GT.DISTM) 

IF ( ISLP1 . LT . ISLPM . AND . 1SLP1 . GT , 800) 

DISTM=DISTS1 
DIST0-DISTS1 
IF( IPE . EO. 1008 ) 

IPED*INT(2 .27>ABS( SIN( ( TIME ( K . I ) » 1  BO  )/( 720'C1 ) ) ) ) 

IPE'IPE-IPED 

END  IF 

IF(TIME(K. I ) .GT.SCMIN( JJ)) 

ISCSLP«SCSLP(OJ) 

ELSE 

CC«FLOAT(SCMIN(vJ)-TIME(K.I)  ) 

DD»FLOAT(SCMIN( J)-SCMIN( J-1 ) ) 

RAT-CC/DD 

ISCSLP«SCSLP( J)-INT(RAT*(SCSlP< J)-SCSlP( J-1 ) ) ) 

END  IF 

IF( LAT(K. I  I  lE  30  ) 
ipn=id:7-:n-(ccsi..at(k.:)'3  /cii-it  p 
ELSE 

IPN« 1027*1NT( COS( LAT(K. I )»3  /C 1 ) ’ 1 7  ) 

END  IF 

IF(IPN-IBE  L'^  3)IPN»IPE-*3 
IF(IPN-ISlP'  l’  3 )IPN»ISLP1*3 
EE«FlOAT( IPN-ISCSLP  ) 

FF. FLOAT! ISLP1-ISCSLP) 

IF( IPE  GT  ISCSLP  AND. ISLP1  GT  ISCSLP) 
AA-(DISTSi-*BB )•( ALOGIEE/FF 1 ' 

GG*flOAT( IPE-ISCSLP ) 

RS'IAA  (A4.0G(EE  GG)))*-)'  'BE)*C3 

IF(ISlf’'  GE  IPE  pIS^pm^ISlp- 

END  IF 

END  IF 

END  IF 

12  CONTINUE 

IF(DISTM.NE  DIST0)RS*DISTM-C3 
IF(DISTP»C3.GT  RS)RS-DISTP*C3 
WRITE! 19. 198)NUM(K).N(K).RS.IPE.IPN.Ni 

98  FORMAT!  IX . 17 . IX . 12. IX . F5  0 . 1 X . 14 . 1 X . 1 4 . 1 X . 1 2  ) 

KLL-0 

KHH«0 

13  CONTINUE 
GO  TO  98 

99  STOP 


107 


END 

MINUTE**»»*»»»»»*»»» 

THIS  SUBPROGRAM  CONVERTS  CODED  TIME  TO  A  FORMAT  THAT  IS  A  CONTIN¬ 
UOUS  MINUTE  STARTING  AT  OOOOZ  ON  THE  FIRST  DAY  OF  DATA  ON  A  PARTI¬ 
CULAR  STORM.  THIS  ASSISTS  IN  COMPUTATIONS  REQUIRING  TIME  OIFFEREN- 
CES  (SYSTEM  VELOCITIES.  ETC.). 

VARIABLE  LIST*** 

MINI  -  MINUTE  REMAINDER 

SUBROUTINE  MINUTE ( ITIM . MIN ) 

MIN*0 

MIN1*ITIM 

100  MIN1«MIN1-10000 

IF(MIN1 -2400)200.600. 600 
600  MIN«MIN+1440 
GO  TO  100 

200  IF  (MINI -100)300. 700. 700 
700  MIN=MIN*60 

MIN1*MIN1-100 
GO  TO  200 
300  MIN=M1N+MIN1 
RETURN 
END 

VEL . * . 

’•HIS  SUBPROGRAM  TAKES  THE  CHANGES  IN  LATITUDE  AND  LONGITUDE  BE¬ 
TWEEN  CONSECUTIVE  FIXES  AND  COMPUTES  THE  RESULTING  SYSTEM  VELOCITY. 

VARIABLE  LIST*** 

DIFLAT  -  CHANGE  IN  LATITUDE 
DIFLON  -  CHANGE  IN  LONGITUDE 
ITIDIF  -  TIME  DIFFERENCE 
HLAT1  -  dummy  argument 
HL0N1  -  dummy  ARGUMENT 
MIND  -  DUMMY  argument 
VELLAT  -  DUMMY  ARGUMENT 
VELLON  -  dummy  argument 

SUBRDUTINE  VEL( TlAT .TLON.MIN.HLATI .HLONI .MIN2 .VELLAT . VELLDN ) 

DIFLAT=TLAT-HLAT 1 

DIFlON*TLON-MLON1 

ITIDIF*MIN-MIN2 

HLAT1*TLAT 

HL0N1*TL0N 

MIN2*MIN 

VELLAT  =  DIFlAT, 'ITIDIF 

VELLON^DIFlON.'ITIDIF 

RETURN 

END 

EXCH**** . 

THIS  SUBPROGRAM  EXCHANGES  FIX  VALUES  SO  THAT  THEY  CAN  BE  USED  AS 
PAST  POSITIONS  AND  TIMES  WHEN  THE  NEXT  FIX  IS  EVALUATED. 

SUBROUTINE  EXCHG(  TLAT , SLAT  1 . TLON . SL0N1 . MIN . MIN2 ) 

SLAT1*TLAT 

SL0N1*TL0N 

MIN2*MIN 

RETURN 


COMP 


THIS  SUBPROGRAM  SEPARATES  THE  OBSERVED  WIND  INTO  U  AND  V  COMPO¬ 
NENTS.  THE  SYSTEM  MOTION  (U.V)  COMPONENTS  ARE  THEN  SUBTRACTED  TO 
ADJUST  THE  WIND  FOR  SYSTEM  MOTION. 

VARIABLE  LIST*** 

FLVEL  -  DUMMY  ARGUMENT 
FUVEL  -  DUMMY  ARGUMENT 
FVVEL  -  DUMMY  ARGUMENT 
IFLANG  -  DUMMY  ARGUMENT 

SUBROUTINE  COMP ( FLVEL . I FLAN6.C1 . SCKLAT . SCKLON . FUVEL . FVVEL) 

IF( IFLANG. EO. 900) 

FUVEL-0. 

FVVEL-0. 

ELSE 

FUVEL*FLVEL*( -SIN{ IF  LANG/C  1 ))-SCKLON 
FVVEL*FLVEL*(-COS( IFLANG/CI ))-SCKLAT 
FLVEL*SORT( FUVEL ••2+FVVEL**2) 

IF( FUVEL. EO.O. . AND . FVVEL . EQ . 0 . ) 

IFLANG*900 

ELSE 

I F  LANG* INT ( ATAN2 ( FUVEL . FVVEL ) 'C 1 -  180 .  ) 

END  IF 
END  IF 

IF ( IFLANG. LE .0)IFLANG*360+IFLANG 

RETURN 

END 

MID . . 

THIS  SUBPROGRAM  CHECKS  INTERMEDIATE  WIND  INFORMATION  THAT  DOES  NOT 
HAVE  PRESSURE  DATA  FOR  INTERMEDIATE  DISTANCE  BETWEEN  THOSE  WITH 
PRESSURE  DATA. 

SUBROUTINE  MID(DISTL1 ,DISTL2 .DISTL3 . I FLG . FUVEL2 . FVVEL2 . A  ) 

IFLG*0 

IF(DISTL2.LT,0ISTL1 .AND.DISTL2.lt. DI5TL3)G0  TO  400 
IF(DISTL2.GT.DISTL1 .AN0.DISTL2.GT.DISTL3)G0  TO  400 
GO  TO  500 
400  A*A*1. 

FUVEL2*0. 

FVVEL2-0. 

IFLG* 1 
500  RETURN 
END 

MAXI* . . . 

THIS  SUBPROGRAM  DETERMINES  THE  MAXIMUM  WIND  BETWEEN  OBSERVATIONS 
WITH  PRESSURE  DATA  AFTER  TRANSPOSITION  AND  WITH  SYSTEM  MOTION 
EFFECTS  ELIMINATED. 

SUBROUTINE  MAXI ( FUVEL 1 . FUVEL2 . FUVEL3 . FVVEL 1 . FVVEL2 , FVVEL3 . FLMAX ) 
FLVEL1*S0RT( FUVEL 1**2*FVVEL 1**2) 

FLVEL2*S0RT(FUVEL2**2+FVVEL2**2) 

FLVEL3*S0RT(FUVEL3**2+FVVEL3**2) 

FLMAX-AMAX 1 ( FLVEL 1 . FLVEL2 , FLVEL3) 

RETURN 

END 

BEE****** . . 

THIS  SUBPROGRAM  COMPUTES  THE  "BB"  SCALING  PARAMETER  USED  IN  THE 


ANALYTIC  MODEL  OF  RADIAL  SLP  PROFILES  OF  TROPICAL  CYCLONES.  IT 
MAKES  USE  OF  THE  INNERMOST  (CORE  REGION)  SLP  GRADIENT. 

VARIABLE  LIST*** 

P  -  TEMPORARY  VALUE  OF  PGRADS 

SUBROUTINE  BEE ( PGRADS. BB) 

P-PGRADS 

IF(P.GT.3.5)P-3.5 
BB*(((P*2. )/3.5)+.5) 

RETURN 

END 
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